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ABSTRACT 

This paper describes the Hubble Space Telescope imaging data products and data reduction proce- 
dures for the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS). This survey 
is designed to document the evolution of galaxies and black holes at z ^ 1.5 — 8, and to study Type la 
SNc beyond z > 1.5. Fivc^ premier inulti-wavclcngth sky regions arc selected, each "with extensive 
multi"wavelcngth observations. The primary CANDELS data consist of imaging obtained in the Wide 
Field Camera 3 / infrared channel (WFC3/IR) and UVIS channel, along with the Advanced Camera 
for Surveys (ACS). The CANDELS/Deep survey covers 125 square arcminutes within GOODS-N 
and GOODS-S, while the remainder consists of the CANDELS/Wide survey, achieving a total of 
~800 square arcminutes across GOODS and three additional fields (EGS, COSMOS, and UDS). We 
summarize the observational aspects of the survey as motivated by the scientific goals and present 
a detailed description of the data reduction procedures and products from the survey. Our data re- 
duction methods utilize the most up to date calibration files and image combination procedures. We 
have paid special attention to correcting a range of instrumental effects, including CTE degradation 
for ACS, removal of electronic bias-striping present in ACS data after SM4, and persistence effects 
and other artifacts in WFC3/IR. For each field, we release mosaics for individual epochs and eventual 
mosaics containing data from all epochs combined, to facilitate photometric variability studies and 
the deepest possible photometry. A more detailed overview of the science goals and observational 
design of the survey are presented in a companion paper. 
Subject headings: Cosmology: observations — Galaxies: high-redshift — 
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1. INTRODUCTION 

In this paper we describe the Hubble imaging and 
mosaic data products from the Cosmic Assembly Near- 
infrared Deep Extragalactic Legacy Survey (CANDELS), 
a 902-orbit Multi-Cycle Treasury (MCT) program aimed 
at documenting the evolution of galaxies and black holes 
from 1.5 to 8, characterizing Type la supernovae 
(SNe la) beyond z > 1.5 to better constrain the nature 
of dark energy, and probing galaxy evolution into the 
epoch of reionization. The CANDELS program uses the 
Wide Field Camera 3 / Infra-Red channel (WFC3/IR) as 
its prime instrument, as well as the WFC3/UVIS chan- 
nel, and obtains parallel observations with the Advanced 
Camera for Surveys (ACS). It is executing across Cycles 
18, 19 and 20, and resulted from the combination of two 
approved proposals that were submitted in response to 
the special Hubble MCT call for proposals 2009, which 
provided for large programs to address unique and broad 
science themes that could not be accommodated within 
the standard annual time- alio cat ion process. 

The structure of the survey includes essential elements 
of two MCT programs that were submitted separately: 
one, led by Ferguson, involved studying the full area of 
the GOOD S-North and GOODS-South fields ( |Giavahsco 
et al.|2004 ) to a uniform depth, including also U V imag- 
ing, and carrying out an extensive search for high-z 
SNe la; the other program, led by Faber, aimed at study- 
ing half the GOODS-North and GOODS-South areas to 
a greater depth, together with wider /shallower imagin; 



a greater deptn, togetner witn wider /snallower irnagmg 
of the Extended Gro th Strip (ECS: |Davis_ei_al, ] [ 2007|, 
Newman etaLj l gOllI in prep) , COSMU5~i |^coviile et al' 
2007l |Koekemoer et al.| | 2007[), and the UKlDSS Ultra- 
Deep Survey (UDS: Lawrence et ai.|2U07 Cirasuolo et al. 
2007 1 , while also permitting a search tor SNe ia. 



2007 1 , while also permitting a search tor SNe ia 

The combined CANDELS program obtains observa- 
tions across all five fields, as well as including the SNe la 
follow-up program, the UV imaging, and the multi-tier 
Deep-fWide observing strategy. Most of the observa- 
tions use WFC3/IR as prime and ACS/WFC in parallel, 
and substituting UV imaging with WFC3/UVIS for parts 
of the orbit in Hubble 's continuous viewing zone (CVZ) 
that are too bright for WFC3/IR imaging. More de- 
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tailed information is presen ted at o ur CANDELS website 
http://candels.ucolick.org and in Grogin et al. (20111 
which provides a more detailed overview of the science 
goals and observing strategy. The outline of this paper 
is as follows: a brief outline of the major science goals is 
given in f2]to place the data products in context; a de- 
scription of the fields is provided in f|3j the observations 
are described in S|4j the data products are presented in 
^ and we conclude in Sj6j 

2. SCIENCE GOALS 

We summarize here the CANDELS science goals in 
the context of h ow they relate to the Hubble data prod- 
ucts. We refer to Grogin et al. ( 2011[ ) for a more detailed 
description of the science goalsj which include stud- 
ies of galaxies in the reionization era ("cosmic dawn"), 
the growth and morphological transformation of galaxies 
during the era of peak star-formation and AGN activ- 
ity ("cosmic high noon"), and measurements of high-z 
SNe la to constrain dark energy and measure supernova 
rate evolution. 



;ions (|Fan et al.| |2006 l 
jpergel et aLpOoVi rm- 
lium was reionized be- 



Cosmic Da wn. QSO observ at 
and WMAP ( |Page et al.|[2007{ [Spe ^ 

dicate that the intergalactic medium was reiomzei 
tween = 0.5 and 1 Gyr. The IGM was seeded with 
metals to Z ^ A x 10""* within the first billion years, 
and the energy released by the stars that produ ced these 
meta, l s appears sufficient to rei onize the IGM (Songaila 



2001 Ryan- Weber et al. 2006). The bright end ot the 



UV luminosity function (U VLF) of star-forming galaxies 
evolves rapidly at 4<z<7 (e.g. , 
Bouwens et aL][2007l [20081 12010| , 

ixies appears msuiiicient 



Dickinson et al.||2004 
but the U V tlux from 
these galaxies appears insufficient to explain reionization 
without extrapolations — an important puzzle to solve. 
Current z ~ 8 luminosity functions are based on only a 
handful of objects, most with Luv ^ L* . The CAN- 
DELS data in the ZS50Y105J125H1Q0 filters can provide 
measurements of the bright end of the UVLF at z 7 — 8 
and also permit robust LBG color selection at (z) = 5.8, 
6.6, and 8.0, for L* for z = 7 (Ji25w27) and 1.5L* for 
z = 8 (i?i6o~27), as well as fainter LBGs at higher z. 
This can constrain extinction via UV spectral slopes and 
improve measurements of their evolution. 

Furthermore, the evolution of faint AGN at 2^6 — 7 
can be directly probed by cross-correlating the drop-out 
samples w ith the d eep X-ray data in these fields (e.g.. 
Fan et al.ir2 003 , Ko ekemoer eTar][2004| |Aird et al.||20q8] 



Brusa et al.,j2009|) for which the depth in the near-IR 
data is cruciaT Moreover, photometric redshifts and 
non-LBG color criteria can help reveal whether there 
are non-star-forrning galaxies lurking at these redshifts 

/n\/r„u„„i, „„ „j- „i I'onnci nir;i.i;„j „4- „i lonnol l/^i,„„. „j- „i 



2008 Chary et al. 



j^obasher et al."2005' '' yiklind et al 

[2D07j Dunlop et al. 2007^), thereby driving the photomet 
ric requirements of the CANDELS data products. The 
CANDELS data also allow fiuctuations in the extragalac- 
tic background light (EBL) to be probed, potentially con- 
straining the proper t ies of the first generati ons of stars 



(jCooray et al. 2004j fFernandez et al. 2010|), as well as 
cf 



enabling the use ot clustering st atistics to constrain the 



properties of dark matter halos (jConroy et al.|2006 
et al.|2006[[2009[ ), which drives t 



Lee 



le need to produce con- 



tiguous mosaics across each of the CANDELS fields. 
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Cosmic High Noon. At z^2, star formation and nu- 
clear activity within galaxies are at their peaks while 
the morphological differentiation of galaxies is well un- 
der way. A key question is what drives stellar mass 
buildup, bulge growth, and the emergence of passive el- 
lipticals. To resolve this requires accurate mass func- 
tion measurements well below M*, achieved by robust 
SED fitting at rest-frame optical wavelengths, where the 
4000 A and Balmer breaks constrain accurate photo-z 
measurements, stellar population ages, and M/L ratios. 
The CANDELS data are designed to match the photo- 
metric depths of the existing Spitzer IRAC and HST ACS 
data, accurately determining the mass function of quies- 
cent galaxies for M > lO^Af0 at z«2 and their contri- 
bution to the global mass density. Another key question 
in galaxy growth is the relative importance of mergers 
and structural instabilities. The discovery of large, ro- 
tating, clumpy, gas-rich disks at z '--^ 2 implies that disk 
instabilities may dri ve bulge formation more rapidly than 



previously thought (Forster Schreiber et al. 2006 2009 



Genzel et a l. 2006, 2008]) 



i'he star formation rate also correlates strongly with 



lation declines steadily below z ~ 


2^Daddi et al.||2007 


Elbaz et al. 


2007 


Noeske et al. 


2007). The empirica 



evidence is now reminiscent of the "clump-merging" see 
nario for the growth of bulges in gas-rich disk s from nu- 



merical simulations (e.g., [Noguchi 1999; Immeli et al 



2004] 



Bournaud et 



TElmcgrcen et ai. ' 2008 ) 



CANDeLS can provide a census of clumps within galax- 
ies along with their sizes and masses, and the result- 
ing estimates of bulge formation rates can be compared 
to the timescale of clump migration driven by dynami- 
cal friction. Through comparison with the deep X-ray 
catalogs, CANDELS also provides detailed morphologi- 
cal information on AGN in this redshift range, tracking 
the connection between galaxy mergers and black hole 
growth. Finally, the evolution of galaxies with very low 
specific SFRs (passive galaxies) is also of interest, with 
sources ha ving been found out to at lea st z = 2.5 (e.g., 
Daddi et al. 2005; Trujillo et al. 2006; [Cimatti et aLl 



2008[ [van Dokkum et al..,2008^ . CANDELS enables the 
luminosity function of large numbers of these sources to 
be directly constrained, in addition to probing their mor- 
phogical structure to faint limits in the ncar-IR with bet- 
ter resolution than previous studies. 

Deep Ultraviolet Observations. An important fea- 
ture of CANDELS is the fact that the GOODS-North 
field is in the HST CVZ, thus we use the bright day-side 
of the orbit to observe with WFC3/UV1S in the ultravio- 
let (F275W and F336W). This enables measurements of 
the Lyman-continuum (LyC) escape fraction from galax- 
ies at z ^ 2.5, identification of ~350 Lyman-break galax- 
ies at z '-^ 2, and measurements of the star-formation rate 
in low-lu minosity dw arfs whi ch may just be "turn ing on" 
at z~ 1 ([Babul & Recs 1992T| Bullock et al.|2000 ]. There 
are ~40-50 UV-luminous LBGs (Luv > 0.25L*) in 
this field at 2.38<z<2.55 (half with spec-z), which is 
the optimal redshift for constraints with the F275W fil- 
ter, many of w hich may be bright enough to detect i f 



/osc>0.5 (e.g., Shapley et al. 2006 Iwata et al. 



Importantly, these galaxies are at redshitts that al 



2om. 

ow Ha 



measurements for an independent measure of the ioniz- 
ing continuum. The resolved LyC distributions provide 
tests of different mec hanisms for high fesc includ ing SNe 



winds ( jClarke fc Oe y 2002, Fujit a et a l.||2002[), minor 
galaxy interactions ( jGnedin et al. 20118), ancTemission 
from globular cluster formation ( ,Ricotti^^2002) . 

Supernova Cosmology. While SNe la at z< 1 — 1.5 
have already provided startling evidence of dark energy 



(iRiess et ai.||1998'; Pcrlmutter et al.||1999l iRiess et al. 



2004 , ' 2007p , CANDELS now provides a direct probe of 
1.5 < z < 2.5 to test the n ature of SNe la progenitors and 
their possible evolution ( Ruiz-Lapuente fc Ca nal''1998[ 
Mannucci et al.||2005l [kobayashi fc JN omoto||2b 09 Greg] 



gio et al.||2008p 7"wEich can be tested withTANDEL* 
since the predicted rates diverge significantly at z > 1.5. 
In addition, CANDELS yields SNe la at z < 1.5 which re- 
main crucial tracers of the evolution in the dark energy 
equation of state w. Increasing the samples at 0.7 < z < 1 
in the IR reduces the uncertainties in host extinction, 
thereby testing whether dust is a factor in the declining 
high-z SNe la rate. CANDELS also includes foUowup 
WFC3 or ACS grism observations to determine the su- 
pernovae type and redshift, and rest-frame B and V light 
curves for each SN (thus limiting the K-correction errors 
to below the random distance error of the SNe). We note 
also that related follow-up programs to CANDELS (in- 
cluding the SNe la search) obtain grism data on these 
fields, but these are separate programs from the CAN- 
DELS imaging survey and are not discussed in further 
detail here. 



3. THE CANDELS FIELDS 

Here we summarize th e general properties o f the CAN- 
DELS fields and refer to Grogin et al. (2011 1 for a more 
detailed description. The CANDELS survey consists 
of a two-tier Deep-fWide survey designed to address 
the science goals discussed in Sp] as well as providing a 
legacy dataset on these fields. The CANDELS Deep por- 
tion covers ~ 125 square arcminutes to ^ 10-orbit depth 
within the GOODS-N orth and GOODS-Sou th Fields 
("Giavahsco et al.|2004[) including the E- CDFS (|Rix et al. 



2004) as well as the WFC3 ERS2 field (| Windhorst et al. 
(2011) 
tal of 



The full area of the CANDELS survey covers a to 
^ 800 square arcminutes, where the additional area 
includes the shallower Wide portion to ~ 2-orbit depth 
around the Deep portions of GOODS, together with sub- 
sections of three a dditional fields, n amely the Extended 
Groth Strip (EGS:' Davis et al.|2007| Newma n et al. 2011, 
in prep.), COSMOS (IScoville et al.||2007| |Koekemoer 



let al. 2007 ),and the UKIDSS Ultra-DeeTS urvey (UPS: 
|Lawrence et al.|2007 Cirasuolo et al.|2 d07|. When com- 



bined with the existing Ultra Deep I' ield (UDF: 



Beckwith 



et al |2006| within GOODS-South, CANDELS provides 
a unifying survey at three principal exposure time depths 
with approximately an order of magnitude difference be- 
tween each. Another unifying aspect of the survey is that 
all five CANDELS fields are the targets of the Spitzer 
Extended Deep Survey (Fazio et al., in prep.) which 
is a 2108 hour program, together with a more recently 
approved 1200 hour follow-on program, covering each of 
these regions with SpitzerlRAC 3.6 /Ltm and 4.5 /zm imag- 
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Table 1 

New HST WFC3 and ACS Filter Coverage in the CANDELS Fields 



Field 


Prog. 


Cycles 


WFC3/IR 




WFC3/UVIS 




ACS/WFC 






ID 




Filters 


#Exp 


Filters 


#Exp 


Filters 


#Exp 


UDS (Wide) 


12064 


18 


F125W F160W 


352 


F350LP 


24 


F606W F814W 


376 


GOODS-S Deep 


12060,1,2 


18+19 


F125W F160W 


696 


F350LP 


174 


F606W F814W F850LP 


870 


GOODS-S Wide 


12060,1 


18+19 


F125W F160W F105W 


134 


F350LP 





F606W F814W F850LP 


134 


EGS (Wide) 


12063 


18+20 


F125W F160W 


360 


F350LP 


18 


F606W F814W 


378 


GOODS-N Deep 


12442,3,4,5 


19+20 


F125W F160W F105W 


660 


F350LP 


360 


F606W F814W F850LP 


780 


GOODS-N Wide 


12442,3,4 


19+20 


F125W F160W F105W 


240 


F350LP F275W F336W 


40 


F606W F814W F850LP 


240 


COSMOS (Wide) 


12440 


19 


F125W F160W 


352 


F350LP 


24 


F606W F814W 


376 


Total 








2794 




640 




3154 



ing to a total depth of 12 hours per pointing. 

Eaeh of the five CANDELS fields has a wealth of ad- 
ditional imaging and spectroscopic ancillary data from 
X-rays to radio wavelengths, described in the aforemen- 
tioned papers and their references. For the present work, 
we note in particular that all of them, except the UDS, 
have pre-existing HST data covering the field. In addi- 
tion, all five fields have extensive pre-existing catalogs 
that can serve as astrometric and photometric reference 
standards as well as being combined with the catalogs 
from new HST data to obtain derived measurements of 
source properties including photometric redshifts, stellar 
masses and star formation histories. 

4. OBSERVATIONS 

We summarize here the general lay out of the observa- 
tions, referring to Grogin et al. (2011 ) for a more detailed 
description. The HST observations in the CANDELS 
fields can be summarized as consisting of three comple- 
mentary sets of imaging data: WFC3/IR, WFC3/UVIS 
and ACS/WFC imaging exposuresj^p] In all cases, the 
WFC3 observations are taken as the prime observations, 
while the ACS data are obtained in parallel. The filter 
breakdown is shown in Table [T] for each of the different 
fields. We also discuss two generally different sets of ob- 
servations, namely the GOODS fields, which contain the 
CANDELS-Deep pointings (together with a CANDELS- 
Wide "flanking field" in GOODS-South), and the other 
three fields (COSMOS, EGS and UDS) which only con- 
sist of the CANDELS- Wide component and thus have a 
somewhat different structure. 

4.1. Filters and Exposure Times 

We first describe the filter choice and exposure struc- 
ture of the three Wide fields (COSMOS, EGS and UDS), 
each of which is covered using a mosaic grid of tiles and 
repeated over two epochs. During each epoch. Each tile 
is observed for one orbit (~ 2000 s), divided into two ex- 
posures in F125W (at a depth of ~l/3 orbit) and two 
exposures in F160W (at a depth of ~ 2/3 orbit), together 
with parallel exposures using ACS/WFC in F606W and 
F814W. However, some WFC3/IR tiles near one end of 
the Wide mosaics are not covered by ACS parallels so 
a short 434 s WFC3/UVIS F350LP exposure is inserted, 
creating a total of five WFC3 exposures per orbit for 
these tiles, and five ACS/WFC exposures in parallel. 

The GOODS-North and GOODS-South fields contain 
the Deep portions of the CANDELS survey, with a total 



For current details on ACS see http://www.stsci.edu/hst/acs 
For current details on WFC3 see http: / /www. stsci.edu/hst/wfc3 



depth of at least 4 orbits in both WFC3/IR F125W and 
F160W and 3 orbits in F105W, spread across 10 epochs. 
Each single-orbit pointing, for each epoch, contains four 
WFC3/IR exposures (two F125W and two F160W), and 
one WFC3/UVIS (F350LP). In parallel, we also obtain 
five ACS/WFC exposures, where the primary require- 
ment is to obtain at least 32,000 s depth in F814W. 
Once this requirement is met, the next ACS/WFC pri- 
orities are - 2500 s of F850LP, followed by - 5000 s in 
F606W, and finally any remaining depth is placed back 
into F814W. Since the GOODS-North field is in the CVZ, 
some portions of the orbit are too bright for observations 
using WFC3/IR so WFC3/UVIS exposures are substi- 
tuted, using the F275W and F336W filters. In these cases 
the ACS parallels retain their structure as described for 
the remainder of the Deep observations. 

Finally, the GOODS-South Deep portion has a wider 
"flanking field", similar to the Wide fields in its filter 
choice and exposure time requirements. This is divided 
into two epochs, achieving a depth of ~l/3 orbit in 
F125W and ^2/3 orbit in F160W per epoch, also us- 
ing short WFC3/UVIS F350LP exposures where nec- 
essary, and with ACS/WFC F606W and F814W expo- 
sures in parallel. For the ACS/WFC exposures falling 
outside the Deep area, the exposure time requirements 
are firstly to obtain ^2500 s in F814W, followed by 
~ 2500 s in F850LP, and distribute any remaining time 
according to the same priorities as for the Wide mosaics. 
The GOODS-North "flanking field" is similar in design, 
with a slightly different layout. Note that, since the 
ACS/WFC field of view is larger than WFC3/IR, the 
ACS parallel pointings overlap considerably and their ef- 
fective exposure time on overlapping pointings can ex- 
ceed the nominal 2 orbit observing time, due to the dense 
packing of pointings for contiguous WFC3/IR coverage. 

In Table [2] we list the current zeropoints (in the AB 
magnitude system) for all the filters used in the CAN- 
DELS survey observations. We note that these are sub- 
ject to change and we provide links to the instrument 
webpages where the most up-to-date zeropoints can be 
obtained in future. The primary uncertainties associated 
with these are related to the spectral characteristics of 
the standard stars that are used by staff at STScI to 
carry out the photometric zeropoint calibrations, as well 
as changes in the instrument and filters over time. Gen- 
erally these are accurate t o be tter than ~ 1 — 2%, and 
we present later in Section |5.9| a quantitative validation 
of this level of accuracy using the photometry from our 
CANDELS data, compared with photometric data from 
ground-based imaging. 
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Table 2 

CANDELS ACS and WFC3 Filter Zeropoints^' 



Instrument/Camera Filter Zeropoint (ABmag) 



ACS/WFC 


F606W 


26.49 


ACS/WFC 


F814W 


25.94 


ACS/WFC 


F850LP 


24.84 


WFC3/UVIS 


F275W 


24.14 


WFC3/UVIS 


F336W 


24.64 


WFC3/UVIS 


F350LP 


26.94 


WFC3/IR 


F105W 


26.27 


WFC3/IR 


F125W 


26.25 


WFC3/IR 


F160W 


25.96 



For current filter zeropoint information, please see 
the following wcbpages: 

Ihttp: / / www.stsci.edu/hst / acs/analysis/zeropoin ts ] 
[|Sttp:// www. stsci.edu/hst/wtc3/phot_zp_lbn 



4.2. Mosaic Layout Design 

For each of the five CANDELS fields, the goal is to 
cover a contiguous area with WFC3/IR (thus, the larger 
ACS/WFC parallel exposures overlap somewhat to cre- 
ate deeper pointings), and to overlap as much as possible 
the existing relevant ancillary datasets. Here we summa- 
rize the specific considerations for each of the fields and 
how they impact the overall design of the mosaic obser- 
vations. 

For the three Wide fields (COSMOS, ECS and UDS), 
the layout consists of a rectangular region, which for 
COSMOS and UDS comprises a grid of 4x11 tiles 
(^8f6 X 23.'8), at spacing intervals designed to allow for 
maximal contiguous coverage in WFC3/IR without in- 
troducing gaps between tiles as a result of pointing errors. 
The exposures are all oriented so that the ACS/WFC 
parallels are offset along the long axis of the mosaic, 
thereby producing a similar-sized mosaic overlapping 
the bulk of the WFC3/IR mosaic, except at its ends 
where some tiles are covered only by WFC3/IR or by 
ACS/WFC, but not both. For the ECS field the mosaic 
is instead 3x15 tiles (~ 6^5 x 32'.5), to optimize coverage 
with ancillary data. 

For the GOODS-North and GOODS-South Deep re- 
gions, the layout consists of a smaller rectangular grid 
of 3x5 tiles (~6;5xl0;8). In GOODS-South the 
WFC3/IR pointings are placed adjacent to the exist- 
ing WFC3/IR ERS2 observations (Windhorst etal 2011). 
Since the field is observed over 10 epochs, the orienta- 
tion of the tiles rotates by ^ 45 — 50° from one epoch to 
the next. This also changes the coverage of the parallel 
ACS/WFC observations, creating a net effect of a larger 
area covered by ACS/WFC to shallower depth, surround- 
ing the deep central WFC3/IR data, which therefore has 
a slight deficit of ACS coverage. 

Finally, the shallower "flanking field" region in 
GOODS-South covers ~ 2 x 4 tiles, divided into two 
epochs, and using a similar filter choice and exposure 
time strategy as the other Wide fields. In Goods-North 
the layout is similar, with differences due to the field 
geometry. The pointings are oriented such that the 
ACS/WFC parallels land mostly on the GOODS-South 
Deep region. The GOODS-South "flanking field" region 
is also covered with F105W to 1 orbit depth which pro- 
vides additional parallel ACS/WFC data for the central 
region. 



Table 3 

CANDELS WFC3/IR Sub-pixel Dither Pattern 



Position 


Filter 


X offset (pix) 


Y offset (pix) 


1 


F160W 


-2.5 


-2.5 


2 


F125W 


-4.5 


-1-2.0 


3 


F125W 


-1-2.0 


-1-2.5 


4 


F160W 


-1-4.0 


-2.0 



4.3. Sub-Pixel Dither Pattern 

Each mosaic tile is observed for one orbit during each 
epoch, where the prime WFC3/IR observations consist 
of four exposures. In most cases these four exposures 
consist of two exposures each in F125W and F160W, ex- 
cept for the Y-band visits where all four exposures are 
obtained in F105W. Due to the relatively large pixel scale 
of the WFC3/IR detector (O'.'128/pixel at its central ref- 
erence pixel), we offset the four WFC3/IR exposures in 
each orbit using a 4-point small-scale dither pattern to 
provide half-pixel subsampling of the point spread func- 
tion (PSF) while also mitigating the impact of hot pixels 
and persistence. The strategy of dividing the 4-point 
dither pattern into two pointings with F125W and two 
with F160W was particularly motivated by the super- 
nova science, where we carried out tests to ensure that 
good supernova subtraction could be achieved. Since the 
second epoch on these tiles also contain two exposures 
in F125W and F160W, we ultimately achieve a 4-point 
dither pattern in each filter, on all tiles, once the data 
from all the epochs are combined. 

The dither pattern serves two complementary pur- 
poses: 1) provide non-integer shifts to subsample the 
PSF; 2) add integer components to these shifts in order 
to ensure that hot pixels and possible persistence from 
previous bright sources are moved around sufficiently. 
In particular, persistence is a concern, especially "self- 
persistence" from sources in previous CANDELS expo- 
sures executed as part of the dither pattern, since it typi- 
cally tends to be more extended than a single pixel and is 
diffuse, thereby subtly impacting the photometry if it is 
not mitigated. For compact sources, the expected spatial 
extent has been quantified as a circle ~ 2.5 pixels in diam- 
eter (WFC3/IR), and therefore constrains the minimum 
size of the dither offset. On the other hand, due to the 
geometric distortion of the detector, an offset of a certain 
number of pixels at the center corresponds to a different 
number of pixels near the edge, and for sufficiently large 
shifts the subsampling can vary from half-pixel to inte- 
ger several times between the center and the edge of the 
detector, introducing non-uniform sub-pixel sampling. 

Therefore, the desire to retain uniform half-pixel sam- 
pling across the entire WFC3/IR detector during each 
epoch constrains the dither offsets to be as small as pos- 
sible, which in this case is the minimum size needed to 
avoid issues from persistent regions ~ 2.5 pixels across. 
As a result, larger detector blemishes (such as the 
WFC3/IR "deathstar", a circular region of bad pixels 
^ 50 pixels in diameter, or ^ 6'.'4) are not covered during 
a given epoch, and instead subsequent epochs at differ- 
ent orientations are used to provide coverage. The final 
pattern that satisfies both the persistence requirements 
and the PSF sampling requirements is a 4-point dither 
pattern with parameters shown in TablejSj and presented 
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Figure 1. The CANDELS WFC3/IR sub-pixel dither pattern, 
indicated by the set of four red labels HI, J2, J3, H4, as well as 
the standard default dither pattern (indicated by the four black la- 
bels). The large circles indicate the 2.5-pixel diameter regions that 
are expected to be impacted by persistence from a point source, 
thus their overlap ends to be minimized. The new CANDELS 
4-point dither pattern is chosen to be a few pixels wider, simul- 
tanteously still providing half-pixel subsampling for the WFC3/IR 
observations as well as introducing larger offsets to help mitigate 
the effects of persistence from a previous exposure, by ensuring 
that any persistent pixels from previously observed bright sources 
are moved around to different, non-repeated pixel locations and 
can therefore be excluded from the final image combination. 

graphically in Figure [l] This has the advantage of being 
large enough to avoia small-scale detector defects and 
persistence from compact sources, but small enough for 
the half-pixel subsampling to change by < 0.1 — 0.2 pixels 
across the detector. 

4.4. Detector Characteristics and Read-out Modes 

All the near-infrared HST observations on the CAN- 
DELS fields are obtained using the WFC3/IR detector 
which consists of a 1024x1024 Teledyne HgCdTe array, 
comprising a central 1014x1014 pixel area for imaging 
(covering a region ^ 130" across, with a plate scale of 
O'.'128/pixel at its central reference pixel), surrounded by 
a 5-pixel wide strip of reference pixels along the edges 
which are unexposed to light and are used to track 
changes in bias level that may occur during an expo- 
sure. The exposures are all obtained using the WFC3 
IR-FIX aperture, which samples the full imaging field of 
view of the WFC3/IR detector. All exposures are ob- 
tained in MULTIACCUM mode using either a SPARS50 
or SPARSlOO read-out sequence, enabling the pixels to 
be sampled non-destructively every 50 or 100 seconds 
respectively. These SPARS read samples are repeated 
for each sequence by specifying the NSAMP parame- 
ter, which range from 9—14 for the SPARS50 sequences 
and 6—10 for the SPARSlOO sequences before reading 
out the array, depending on the scheduling constraints 
for each particular field and epoch, leading to exposure 
times in the range 450 — 1000 seconds. In addition, each 
of these MULTIACCUM exposures is preceded by a short 



read 2.9 seconds in duration (the zeroth read), which 
serves as a measure of the bias structure across the array 
at the start of each exposure. 

The optical and UV observations are all obtained us- 
ing the WFC3/UVIS and ACS/WFC cameras, which 
are rather similar in construction, each consisting of two 
thinned, back-illuminated CCDs with a usable area of of 
either 4096x2048 (ACS) or 4096x2051 (WFC3), located 
adjacent to one another to form an imaging area approxi- 
mately square in size, with a small physical gap between 
the detectors. The WFC3/UVIS wavelength response 
is optimized to be UV sensitive, as compared with the 
ACS /WFC which we use for observations at longer wave- 
lengths. Each detector has two amplifiers, one at each 
corner, with half the pixels being read out in parallel by 
each. The ACS/WFC CCDs were manufactured by SITe 
(Scientific Imaging Technologies) and are each physically 
4144x2048 pixels. The first and last 24 columns are read 
out as physical overscan pixels which are unexposed to 
light, while an additional 20 rows of virtual overscan are 
added while clocking out each column during read-out. 
This yields a usable image area of 4096x2048 pixels for 
each of the two chips, provided by the WFC aperture 
which we use for all the ACS observations, with the phys- 
ical gap between them corresponding to about a 50 pixel 
width (~2'.'5), and covering a combined area ~200" in 
extent with a plate scale of O'.'05/pixel at the central ref- 
erence pixel. The WFC3/UVIS CCDs were manufac- 
tured by e2v (formerly Marconi) and are each physi- 
cally 4146x2051 pixels, being read out using the first 
and last 25 columns as physical overscan pixels which 
are unexposed to light, while an additional 30 columns 
of serial virtual overscan are added at the boundary be- 
tween the two amplifier regions on each CCD, together 
with an additional 19 rows of parallel virtual overscan 
added at the end of each column during read-out. Thus 
their usable image area, provided by the UVIS-IR-FIX 
and UVIS-CENTER apertures which we use for all the 
WFC3/UVIS observations, is 4096x2051 pixels for each 
of the two chips, with the physical gap between them 
corresponding to about 30 pixels (~1'.'2), and covering 
a combined area ^ 160" in extent with a plate scale of 
O'.'0396/pixel at the central reference pixel. 

The WFC3/UVIS camera is used for the UV F275W 
and F336W observations in the CVZ for GOODS-North, 
as well as for the F350LP observations for all fields, which 
is a long-pass filter covering the detector response from 
~ 3500 A to its red cutoff around 1 ^m. The latter expo- 
sures are all 434 s in length, with no CRSPLIT, consisting 
of only a single exposure per visit. The ACS exposures 
are aU obtained in the F606W, F814W and F850LP fil- 
ters, depending on which field is being observed, and are 
all obtained as parallel exposures to either the WFC3/IR 
or WFC3/UVIS exposures. As a result, they have a 
range of exposure times between 225 s and 807 s, driven 
by read-out times and overhead considerations associated 
with managing the buffer dumps of both ACS and WFC3 
in parallel, depending on which field is being observed. 
In addition, the ACS/WFC coverage from one pointing 
to the next is also less homogeneous, both in terms of the 
filters used and also its depth, particularly for the three 
CANDELS Wide fields, due to the scheduling and buffer 
time management constraints. 
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5. DATA PROCESSING AND PRODUCTS 

For each of the five CANDELS fields, the final data 
products consist of a set of mosaics for each camera/filter 
combination, processed using the "MosaicDrizzle" 
pipeline (described here), which carries out astrometric 
registration and image combination usi ng MultiPrizzle 
( Koe kemoer et al. 2002), and Drizzle (Fruchter & Hook 
2UU2| . Here we describe the original da!ta, together with 
each of the calibration and pipeline processing steps in- 
cluded in the "MosaicDrizzle" pipeline (outlined in 
Figure [2|, as well as the properties of the final com- 
bined mosaic images. All the processing is done on a 
dedicated cluster of linux machines at STScI, using the 
Pyraf/STSDAS packages. 

For all exposures in each visit, the raw files are pro- 
cessed first through a number of basic calibration steps 
that are needed to produce a final calibrated file for 
each exposure. Since these steps are somewhat differ- 
ent for the WFC3/IR detectors than for the CCD cam- 
eras, we discuss the WFC3/IR data calibration steps sep- 
arately. The WFC3/UVIS and ACS/WFC calibration 
steps, however, are very similar to one another and are 
therefore discussed together, with differences highlighted 
where necessary. 

5.1. WFC3/IR Detector Calibration 

Here we describe the calibration steps carried out to 
convert the WFC3/IR data from raw counts to a set of 
final, flatfielded and flux-calibrated exposures, which are 
subsequently all combined to create the final mosaics. 

5.1.1. Standard Pipeline Calibration with. calwfS 

The raw WFC3/IR exposures, still containing all the 
separate MULTlACCUM read samples, are calibrated 
using the Pyraf/STSDAS task calwf Q This task pop- 
ulates the bad pixel arrays using knownbad pixel tables, 
followed by a bias subtraction for each read, based on 
the mean value of all the pixels in the 5-pixel wide refer- 
ence pixel region. It then carries out a subtraction of the 
zeroth read in order to remove the bias structure across 
the detector, followed by a subtraction of the dark cur- 
rent reference files for either the SPARS50 or SPARSlOO 
read-out sequences, as applicable. This is followed by the 
non-linearity correction and photometric keyword calcu- 
lation, using the current filter throughput tables and de- 
tector quantum efficiency curves. 

Once all the separate MULTlACCUM reads are cali- 
brated as described above, they then move through the 
next steps in calwf 3, namely up-the-ramp slope fitting 
and cosmic ray rejection. For each pixel, this step per- 
forms a linear fit to the accumulating counts that are 
sampled during each MULTlACCUM read, while reject- 
ing outliers from the fit as being due to cosmic rays. 
Currently a threshold of 4 sigma is used for this rejec- 
tion, and flags are populated in the data quality arrays 
corresponding to the read during which the cosmic ray 
occurred. A final count-rate value is then computed for 
each pixel using only the unflagged reads, and is stored as 
the count-rate in the final calibrated exposure, while the 
uncertainty in the slope of counts versus time is stored 

5"' Further documentation for all the PyRAF/STSDAS data re- 
duction software is provided at http:/ /stsdas. stsci.edu/ 



in the error extension of the image. The final steps of 
calwf 3 include multiplicative corrections for the detec- 
tor gain and the flatfield structure across the detector, 
appropriate to each filter. Other factors affecting the 
fluxes include different pixel sizes due to detector distor- 
tion and are still present at this stage; these are corrected 
later when the geometric distortion is removed. 

5.1.2. Additional Corrections for Persistence, Warm Pixels 
and Flat Field Residuals 

In addition to the default calibrations, there are also 
a number of additional corrections carried out to further 
improve the WFC3/IR data, but which are not part of 
the standard calwf 3 pipeline. The CANDELS team has 
implemented these as additional steps in our automated 
image processing pipelines, and they are executed in con- 
junction with the standard calwf 3 calibration steps. 

The first of these additional corrections concerns the 
presence of persistent flux in certain pixels due to bright 
sources having been observed in previous exposures, 
which can be a significant issue for the WFC3/IR de- 
tcctorPj For the first few epochs of this program, darks 
from tne WFC3 calibration program that execute just 
prior to CANDELS visits are used to aid in identifying 
and measuring problematic pixels. Pixels with persis- 
tent flux are then identified in these dark frames if they 
exceeded a count-rate threshold of five sigma above the 
mean, and are flagged in the following science exposures. 
For subsequent orbits during a visit, the CANDELS team 
can determine directly from the preceding CANDELS ex- 
posures which pixels may contain sufficient flux to cause 
persistence; the calibration darks are only used for the 
flrst orbit in a visit, when the previous data may be from 
another program and not necessarily accessible. During 
the initial calibrations for the CANDELS project, pixels 
flagged as having signiflcant persistent flux are excluded 
when performing the mosaic combinations, but work is 
being done to see if the persistent flux can be modeled 
and subtracted instead, which would enable these pixels 
to be used. 

Another correction that is implemented at this stage 
is the identiflcation of additional "warm" pixels in the 
exposures, which might be fluctuating and therefore per- 
haps not present in the calibration reference dark flies 
but only in the images. The calibration darks are used 
to identify these pixels if they exceed a threshold of flve 
sigma above the mean, in which case they are flagged 
in the data quality arrays that are associated with each 
exposure, and are excluded from the flnal image combi- 
nation. 

Finally, the initial WFC3/IR data obtained for the 
CANDELS program were calibrated using an early gen- 
eration of flatfield files in the STScI archive pipeline that 
did not fully correct for all the flatfield features present in 
the data. Specifically, these included the IR "blobs" that 
have appeared in the WFC3/IR channel since launch and 
were not present in the ground flats, as well as a ~ 3 — 4% 
residual large-scale variation in the overall structure of 
the flatfleld flies. We therefore created residual flatfleld 
flies which we applied in our pipelines to correct these 
effects after having run the data through the initial cal- 
ibration. Both of these issues have subsequently been 

http: / / www.stsci.odu/hst /wfc3 / ins_performance/persistence / 
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Collect incoming data 

' Monitor dedicated directory for "incoming" data from OPUS 
■ Organize all exposures for ACS, WFC3/IR,UVIS for each visit 
' Rename exposures according to visit / filter / exposure number 



RAW EXPOSURES: 
ACS, WFC3/UVIS 



Calibrate ACS & WFC3/UVIS: calacs, calwfi 

(run separately on each exposure) 

• Populate bad pixel / data quality arrays 

• Bias subtraction 

• Dark current subtraction 

• Flatfield and gain correction 



RAW EXPOSURES: 
WFC3/IR 



CALIBRATED 
EXPOSURES: 
WFC3/UVIS 



CALIBRATED 
EXPOSURES: 
ACS 



Calibrate WFC3/1R exposures: calwf3 

(run separately on each exposure) 
Populate bad pixel / data quality arrays 
Bias subtraction 
Zeroth-read subtraction 
Dark curretit subtraction 
Non-linearity correction 
Up-the-ramp cosmic ray rejection 
Flatfield and gain correction 



I Bias Striping Correction | 



CALIBRATED EXPOSURES: 
WFC3/IR 



CTE Correction 



Persistence Flagging 



HE 



Amplifier Quadrant 
Correction 



Si 



IE 



HE 



Warm Pixel Flagging 



Cross-Correlation Shifis 

Create single-drizzled image for each exposure in an 
orbit, for each instrument 

Cross-correlate all exposures relative to first exposure 

in the orbit, for each instrument 

Fit the Fourier transform peak to obtain shifts 

Apply shifts to input exposures 

Reneat above for 2nd and 3rd iteration 



Manual Quality Validation 

Inspect combined images for: 

- Satellite trail defects 

- Detector issues (quadrant offsets, 
crosstalk etc) 

Valildate astrometric accuracy: 

- Overlay catalogs 

- Inspect dx,dy plots 



RELATIVE ASTROMETRY 
CORRECTED 



Initial Multidrizzie (for each orbit) 

For each exposure, calculate sky, inverse variance 
For each camera/filter combination: 

- Run Multidrizzie single-drizzle images 

- Create clean median image 

- Transform median back to input exposures 

- Run cosmic ray rejection on each input exposure 

- Re-run Multidrizzie new combined image 

- Run catalog on new combined image 

Refine absolute astrometry for all images in this orbit 
Run final catalog for each camera/filter in this orbit 



ABSOLUTE ASTROMETRY 
CORRECTED 



Incorporate each new orbit into mosaics 

For each new orbit, re-mn Multidrizzie to add into 
existing mosaic for this epoch: 

- "Single-Epoch" mosaics 

Compare previous epochs to new epoch to improve 
cosmic ray rejection in all overlapping exposures 
Combine new epoch into previous epochs to create 
multiple mosaics for each of the 5 fields: 

- "Wide-only" mosaics (uniform shallow depth) 

- "Deep" mosaics (relevant in GOODS-N,S) 

- "Full" mosaics (all exposures) 



Figure 2. General overview of the "MosaicDrizzle" pipeline used for the CANDELS data processing, showing each of the steps that are 
used to process each new dataset as it arrives, including automated calibration and MultiDrizzle combination, as well as data quality 
validation, astrometric registration and mosaic combination. 



corrected by constructing an extensive set of sky flats 
for each of the different fihers, and the resulting cor- 
rected flatfield reference files are now in the STScI archive 
pipeline and are used when the calwf 3 calibration soft- 
ware is run. 

We have verified that the degree of flatness of the fi- 
nal sky backgrounds in the images is within ~ 1 — 2% of 
the mean sky level. In particular, after adoping the im- 
proved WFC3/IR flatficlds that we have discussed here, 
we verified that the photometric repeatability across the 
detector was consistent to the same level of accuracy. 



with no significant deviations found beyond ^ 1 — 2%. 
We present a more detailed comparison in Section |5.9| 
between our WFC3/IR photometry and existing ground- 
based photometry, to illustrate the degree of consistency 
that is achieved. 

The final, calibrated images for each WFC3/IR expo- 
sure are then used in the subsequent steps of astrometric 
alignment and MultiDriz zle mosaic combination, as de- 
scribed in more detail in §5.8[ 

5.2. WFC3/UVIS and ACS/WFC Detector Calibration 
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Figure 3. Example showing the removal of bias stripe and CTE deg radation in the AC S data: (left): original exposure after calibration 
with calacs; (middle): sai ne exposure after removal of the bias stripes ( [Grogin et al.|2010[ |; (right): final dataset after subsequent correction 
for the CTE degradation l |Anderson fc Bedin|20"Tol l. 



In this subsection the WFC3/UVIS and ACS/WFC 
observations are described together since both of these 
cameras have rather similar physical characteristics, with 
differences highlighted where necessary. 

5.2.1. Standard Pipeline Calibration wtth calwf3 and calacs 

Each of the raw ACS/WFC and WFC3/UVIS expo- 
sures are initially calibrated using the Pyraf/STSDAS 
tasks calacs and calwf 3 respectively, for which the stan- 
dard steps are quite similar, although additional steps 
are incorporated as described later in this section. For 
the standard calibration, in each case the CCD data first 
have an overall bias correction applied, which is measured 
from the physical and virtual overscan regions that are 
around the edges of the detector. This is done by carry- 
ing out a linear fit along the overscan rows and columns, 
to capture any potential slope across the rest of the detec- 
tor that might be the result of a bias drift. Subsequently 
a bias reference file image is subtracted from each ex- 
posure to account for the remaining pixel-to-pixel bias 
structure in the detectors. In addition, the dark refer- 
ence file is subtracted from each exposure to take into 
account the dark current structure across the detector 
and help mitigate warm pixels that may be present in the 
images. Finally, the multiplicative gain correction and 
flatfield structure reference files are applied, followed by 
photometric keyword calibration using the current filter 
throughput curves and detector sensitivity information, 
thereby resulting in a set of exposures for WFC3 /UVIS 
and ACS / WFC that are calibrated according to the stan- 
dard pipeline calibration. 

5.2.2. Additional Corrections for Bias Striping, CTE 
Degradation and Amplifier Quadrant Ojfsets 

In addition to these default calibrations, a number of 
other corrections need to be carried out particularly on 
the ACS/WFC images, primarily related to the length 
of time that the detectors have been on orbit, and the 
changes in the detector read-out electronics as a result 
of the new CCD Electronics Box Replacement (CEB-R) 
that was installed during Servicing Mission 4 (SM4) to 
restore the instrument to operation. These additional 
steps, described below, are implemented in our auto- 



mated CANDELS image processing pipeline and are ex- 
ecuted in conjunction with the standard calacs calibra- 
tion steps, with some modifications. 

The first of these additional effects is the removal of 
bias striping noise from the ACS/WFC exposures. This 
effect is introduced by the SIDECAR ASIC electron- 
ics (System Image, Digitizing, Enhancing, Controlling, 
and Retrieving / Application Specific Integrated Circuit, 
manufactured by Teledyne) which was part of the new 
CEB-R installed on ACS during SM4. This circuitry 
exhibits a low-frequency 0.001 — 1 Hz) noise with a 
power spectrum similar to 1/f noise, which is introduced 
into one of the reference voltages for the ACS/WFC 
CCDs and subsequently manifests itself as a bias am- 
plitude variation from one row to the next. However, 
this noise has the property of being relatively uniform 
across the detectors, with its amplitude distribution be- 
ing Gaussian with a = 0.75 e~ (significantly less than the 
~3 — 4e~ readnoise of these CCDs). The relative uni- 
formity of the signal along each row enables this striping 
pattern to be characterized an d largely removed using an 
algorithm ( Grogin et al.||20To| which fits each CCD row 
independently, determinmg the background level using 
an iterative cr -clipping technique to reject non-sky pix- 
els, followed by a hybrid mean and median estimator to 
compensate for lower-level signal from faint sources and 
cosmic ray hits. This algorithm has been implemented 
in our CANDELS image calibration pipeline in conjunc- 
tion with the calibration stage and reduces this noise to 
< 0.3 — 0.4 e~ , thereby largely removing its impact on the 
data. 

The second detector effect addressed in our pipelines is 
the impact of Charge Transfer Efficiency (CTE) degra- 
dation. This effect comes about because of the read-out 
scheme that is implemented for CCDs, whereby charge 
that has been detected by each pixel on the array is first 
transferred down all the remaining pixels in the same col- 
umn, and subsequently across all the remaining columns 
to the amplifier where it is read out. As the charge pack- 
ets are transferred from one pixel to the next, charge 
traps that are present in the pixels due to impurities and 
crystalline defects can capture some of the electrons, re- 
leasing them after a short period of time. This leads to 
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a loss of flux in the original pixel and manifests itself as 
deferred-charge trails along the columns behind bright 
pixels in each exposure, while also producing a net as- 
trometric shift up along in the column for bright sources. 
The effect becomes increasingly severe for pixels furthest 
from the amplifiers, which for these detectors are the 
pixels near the chip gaps. The charge trap population 
is created by cosmic ray bombardment and increases de- 
pending on the length of time for which a detector has 
been on orbit, which is now > 8 years for the ACS/WFC. 
While WFC3/UVIS was installed more recently, indica- 
tions are that its CTE is degrading more rapidly than 
ACS and it may also need to be corrected, though its 
level of CTE degradation is currently small and has not 
yet been characterized as fully as ACS, which we focus 
on here. Figure [3] shows an example of the impact of the 
bias stripe and CTE degradation in the images. 

We implement a correction for this CTE degradation in 
the CANDELS imaging pipelines, making use of a pixel- 
based algorithm that carries out a fit to the deferred- 
charge trail behind each bright pixel and redistributes the 
charge into the pixels in a way that is intended to repre- 
sent the original flux as det ected, before any CTE de gra- 
dati on (as implem ented by Anderson fc Bedin||2010 see 
also Massey 2010). A key point about this algorithm is 
that it is eifectively a deconvolution, by virtue of the fact 
that it restores the charge profiles of pixels along a par- 
ticular column to their original shape, which is sharper 
and more concentrated than the observed profiles which 
have been smeared by the deferred charge trails. As such, 
the pixel-to-pixel noise in the final reconstructed image 
is also somewhat higher than in the original exposure. 
In addition, the algorithm includes the effect of read- 
out noise, which is introduced at the amplifier and does 
not participate in the deferred-charge smearing on the 
detector pixels, by treating this as a separate noise com- 
ponent from the noise on the detector. Tests to date 
have shown that this algorithm correctly reproduces the 
expected noise that would be present in the images if 
no CTE degradation had been present, and in addition 
restores both the photometry and the astrometric ac- 
curacy to levels that are comparable to images without 
CTE degradation. Therefore, this algorithm is included 
in our pipelines in conjunction with the standard detec- 
tor imaging calibration steps. We will continue to review 
the performance of this correction as more data are col- 
lected. 

Finally, the ACS/WFC detector amplifier quadrants 
exhibit additional bias-related offsets between one an- 
other that are not fully corrected during standard cali- 
bration. We implement a routine that fits for the differ- 
ences between these, using an iterative clipping proce- 
dure to eliminate signal from astronomical sources and 
preserve only the background fiux, which is then used to 
remove the residual amplifier quadrant differences and 
place all four quadrants on a uniform background level. 

The resulting calibrated, flatficldcd exposures for all 
the CCD observations are then used in the subsequent 
steps of alignment, cosmic ray rejection and image com- 
bination, as described in the next section. 

5.3. Relative Astrometry and Distortion 

Once all the individual WFC3/IR, WFC3/UVIS and 
ACS/WFC exposures have had their detector-level cali- 



brations applied, they are next passed through the rest 
of the CANDELS calibration and mosaicing processing 
pipelines. The overall design of this pipeline is shown in 
Figure [5] and is described here in more detail. At the 
starting point in this process, the data from all three de- 
tectors, WFC3/IR, WFC3/UVIS and ACS/WFC have all 
been calibrated and are in a very uniform set of formats, 
effectively representing the incoming photons detected on 
each pixel. Therefore, from this stage onwards they are 
processed through very similar steps in the remainder of 
the pipeline that generate the higher-level products, and 
are all discussed together in this section, with occasional 
differences highlighted as appropriate. The first stage of 
this pipeline processes all the exposures in a given visit, 
for each of the different cameras, and addresses the rel- 
ative shifts between exposures in each given single-orbit 
visit. 

5.3.1. Dithering and Pointing Uncertainties 

The four WFC3/IR exposures, together with their cor- 
responding ACS parallels, are all obtained in a 4-point 
dither box pattern aimed at providing half-pixel sub- 
sampling for the WFC3/IR detectors along both axes 
of the pixel, together with a small integer pixel shift to 
ensure that bad pixels and other artifacts (e.g., persis- 
tence) are moved around. Due to the geometrical dis- 
tortion of the detector, shifts that are too large corre- 
spond to a substantially different number of pixels along 
the edge than at the center where the shifts are defined, 
which would cause the pixel subsampling phase to change 
across the detector. Therefore the shifts are kept small 
enough to retain the intended half-pixel subsampling 
across much of the detector with no significant change 
in phase. This means, however, that large artifacts such 
as the WFC3/IR "deathstar", or the WFC3/UVIS and 
ACS/WFC detector chip gaps (~ 2 — 3" across) are not 
covered by the dither pattern in a single orbit, there- 
fore we rely on observations during subsequent epochs to 
cover these missing areas and provide coverage across 
the entire field. While the dithers are small, there 
is nonetheless a spacecraft positioning error associated 
with these small angle maneuvers, which is on the level 
-3-5mas (i.e., - 0.05 - 0.1 pixel for the WFC3/UVIS 
and ACS/WFC CCDs). 

In addition to the slight positioning uncertainties in- 
troduced by the small angle maneuvers, an optical off- 
set is introduced whenever a different filter is inserted 
into the optical path. This offset is currently not in- 
cluded in the geometric and astrometric information for 
the HST instruments, thereby leading to an apparent 
astrometric change which can be quite significant, on 
the order ^ 0.2 — 0.3 pixel depending on the filter and 
instrument. Moreover, during each orbit the spacecraft 
undergoes thermal expansion and contraction ("breath- 
ing") due to changes in solar illumination, which lead 
to changes in the optical path length to the detectors, 
hence resulting in slight scale changes from one exposure 
to the next. This scale change is distinct from that which 
is produced by velocity aberration related to the orbital 
motion of the spacecraft around the earth (which is accu- 
rately known and can be corrected) , and therefore needs 
to be treated separately. Finally, uncertainties in guide 
star reacquisition from one orbit to the next can lead to 
errors in position as well as small rotation uncertainties. 
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Figure 4. Distribution of uncertainties in the measured shifts obtained by cross-correlation, for all the exposures for a given instrument 
during a given orbit, for all the data obtained so far using ACS/WFC (left) and WFC3/IR (right). The measured uncertainties are typically 
~ 0.01 — 0.02 pixel, thus ~0.5 — 1 mas for ACS/WFC and ^1 — 2 mas for WFC3/IR due to its larger pixel size. This technique was used to 
correct for the shifts between exposures as a result of spacecraft pointing uncertainties in executing dither offsets (up to ~ 5 mas) as well 
as the more significant astrometric shifts introduced by filter changes (~ 10 — 25 mas, depending on the instrument/filter) which are not 
accounted for by the distortion information and need to be solved for empirically. 



while a full acquisition of a new guide star has astromet- 
ric uncertainties of ~ C/.'S — 0'.'5 (reflecting the absolute 
astrome tric uncertainties i n the Guide Star Catalog 2, 
GSC-2: [Lasker et al.|[2008l ). 



5.3.2. Distortion Calibration Models 

In addition to the pointing and guidestar uncertainties 
related to the spacecraft, the astrometric accuracy also 
depends on the degree to which the detector geomet- 
ric distortions are calibrated. For the ACS/WFC cam- 
era, the detector distortion is modeled using a 4th order 
polynomial which has been measured separately for each 
filter, including also a time- dependent evol ution of the 
first-order linear skew terms ('Anderson'2007) which is in- 
dependent of filter. In addition, the measured distortion 
on the detector deviates slightly from a polynomial de- 
scription, and is captured in a "distortion residual image" 
which is combined with the polynomials when correcting 
for the distortion. For ACS/WFC, the global accuracy of 
the distortion solution is now within ~ 0.02 — 0.03 pixel 
(Anderson 2007). 

For the WFC3/IR and WFC3/UVIS detectors, the on- 
orbit distortion is less well characterized due to the rel- 
atively short time for which the instrument has been 
in operation. The latest distortion solutions have been 
delivered on 11 October 2010 for both cameras (see 



Kozhurina-Platais et al. 2009 for general reference) and 
contain empirical measurements of the distortion in all 
the filters used in the CANDELS program. In all cases 
the distortion is modeled as a 4th order polynomial, but 
there is not yet any evidence of possible skew term evo- 
lution, and there are also not yet any distortion residual 
images that are to be combined with the polynomials. 
The global accuracy of the present solutions is < 0.1 pixel 
in each camera, as verified also by our testing on the 
exposures. If these solutions improve in future, we will 
include them in subsequent reprocessing of the data. 



5.3.3. Cross- Correlation Shift Determination 

To solve for and remove the residual uncertainties in 
the spacecraft dither offsets between all the exposures 
in each orbit, a cross-correlation prcedure is applied to 
all the exposures for a given instrument in each or- 
bit. The involves taking the first exposure as a refer- 
ence, and cross-correlating all the subsequent exposures 
against that. Filters that are close in wavelength (e.g., 
F606W and F814W for ACS, or F125W and F160W for 
WFC3/IR) can be successfully cross-correlated because 
the strength of the cross-correlation signal from even a 
small sub-pixel shift is coherent across the image and 
is much stronger than any morphological differences in 
galaxies between adjacent filters. Thus, for example, if 
an orbit contains four WFC3/IR exposures (two F125W 
and F160W), then the first exposure would be used as a 
reference and the other three would be cross-correlated 
against it (since they are all obtained with the same in- 
strument, even though the filters are different). Simi- 
larly, if the same orbit contains four parallel ACS expo- 
sures, then the first ACS exposure would be used as a 
reference and the other three would be cross-correlated 
against it; the ACS and WFC3 are treated completely in- 
dependently since the detectors have quite different prop- 
erties. Currently this technique is only applied within a 
single orbit since other astrometric errors such as orien- 
tation differences can be present between different orbits, 
or different visits, and these are solved for in a different 
part of the pipeline. 

The cross-correlation procedure first passes all the ex- 
posures for each instrument (using all the filters) through 
a partial run of MultiDrizzle, up to the point where 
single-drizzled images are produced for all the individual 
exposures. These images have had the instrument dis- 
tortion removed, as well as having had the commanded 
spacecraft offsets applied, and are all aligned on the 
same pixel grid so that astronomical sources should be 
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at the same pixel locations if no residual shifts were 
present. These images are then masked so that all empty 
or faint regions are set to zero, retaining only regions 
around objects that contain sufficient signal for the cross- 
correlation. Bright saturated sources, in particular stars 
with long charge-bleed columns, are also excluded since 
such columns degrade the quality of the cross-correlation 
solutions. The regions containing all the remaining ob- 
jects are also tapered to avoid introducing artificial ring- 
ing in the cross-correlation step. 

Each of these images is then cross-correlated relative 
to the first one, producing a Fourier transform image 
for each exposure. The Fourier transform image con- 
tains a strong peak, offset from its center by an amount 
that corresponds to the residual shift between the two 
exposures. The profile of this peak is fit using a two- 
dimensional fitting routine to determine its location and 
associated uncertainty, which is then directly translated 
into a shift between the two exposures. The uncertainties 
are typically less than a few hundredths of a pixel. 

After a set of first-pass cross-correlation shifts has 
been obtained, these shifts are propagated back to the 
input exposures, in order to do a second-pass run of 
MultiDrizzle which includes a cosmic-ray rejection step 
with the improved shifts. The pixels flagged as cosmic 
rays in each exposure then have their flux replaced with 
pixels from the clean drizzled image, then subsequently 
re-drizzled to a new set of single-drizzled images for each 
separate exposure, which are again masked and used as 
input for the second-pass cross-correlation step. This 
provides a much cleaner cross-correlation signal, since 
much of the noise in the flrst-pass cross-correlation is due 
to cosmic rays. Finally, the shifts from the second-pass 
cross-correlation are propagated into the image headers 
for a third-pass MultiDrizzle, cosmic ray rejection and 
cross-correlation step, and are also stored as the final 
set of relative coordinates for all the exposures within a 
given orbit. At this stage the expected shifts are zero, 
hence any remaining shifts between the second and third 
cross-correlation iterations provide a good diagnostic of 
the residual uncertainties in the shift measurement. 

Figure |4] shows the distribution of the uncertainties 
in the measured shifts, for all CANDELS exposures 
obtained to date in WFC3/IR and ACS/WFC. These 
are derived from the uncertainties on the location of 
the cross-correlation peak, which are obtained during 
the fitting procedure. The uncertainties are typically 
~ 0.01 - 0.02 pixel, thus -0.5 - 1 mas for ACS/WFC and 
~ 1 — 2 mas for WFC3/IR due to its larger pixel size. 
Occasionally, the uncertainties are somewhat higher de- 
pending on the structure of the sources within a given 
exposure, or residual anomalies such as satellite trails or 
other defects that might be present within the sources 
used for the cross-correlation, but in all cases they are 
still less than O'.'Ol. This technique was used to correct 
for the shifts from one exposure to the next in each filter 
as a result of pointing uncertainties in executing dither 
offsets (up to ~ 5 mas) as well as the much larger astro- 
metric shifts introduced by filter changes 10 — 25 mas, 
depending on the instrument/filter) which are not ac- 
counted for by the distortion information and need to be 
solved empirically. A comprehensive program of testing 
has been carried out to validate this routine, and its final 
implementation in the pipeline is able to correct the rel- 



ative shift errors present between the exposures in each 
orbit to a level of accuracy better than a few milliarc- 
seconds, thereby correcting the small errors introduced 
when the spacecraft executes small angle maneuvers for 
dither offsets, as well as correcting the offsets introduced 
by filter changes. 

It should be noted that this level of accuracy refers 
specifically to the overall relative alignment from one 
exposure to the next during an orbit as determined by 
cross-correlation, and does not reflect the astrometric ac- 
curacy of individual sources across the images. These 
are still limited by the accuracy of the distortion models 
along with the absolute astrometric uncertainties, which 



are further discussed in ^5.5 



5.4. Cosmic Ray Rejection 

With the relative shifts within each orbit having been 
corrected, the next step consists of creating a cosmic ray 
mask for all the exposures of a given filter, for each cam- 
era, during a given orbit, by carrying out another run 
of MultiDrizzle, this time with the improved relative 
shifts. The cosmic rays are identified in the driz_cr step 
of MultiDrizzle using a process that first creates a se- 
ries of separately drizzled images, one for each input ex- 
posure, which are subsequently used to create a median 
image using the "minmed" algorithm in MultiDrizzle, 
which enables the minimum to be used instead of the me- 
dian in cases where valid pixels from only two or three 
exposures are present, if one of them exceeds the others 
by > 4cr. The clean median image is then transformed 
back to the distorted detector frame of each input expo- 
sure to carry out cosmic ray rejection using the following 
approach. The input counts in a given pixel in the orig- 
inal exposure, /exp are compared with the counts from 
the median image /med for the same pixel, together with 
the derivative of the median image Ajnod: defined as the 
steepest gradient from that pixel to its surrounding pix- 
els (with all these quantities being in units of electrons). 
A pixel is flagged as a cosmic ray if it exceeds a threshold 
defined as follows: 



l/n- 



/„,cd| > Sl^^^d^ SNRJa. 



read 



B\ 



where S and SNR are adjustable scaling factors and B 
is the background sky value that has been measured for 
the exposure. The inclusion of the gradient term Amcd 
effectively "softens" the cosmic ray rejection in regions 
of relatively steep gradients such as bright cores of ob- 
jects, where the pixel-to-pixel variation can exceed sim- 
ple Poissonian statistics. For the data processing in the 
CANDELS pipelines, this rejection is performed over two 
iterations, with the flrst pass going through all the pixels 
in the image and using 5* = 1.2 and SNR~3.5, followed 
by a second pass in a 1-pixel wide region around each of 
the pixels flagged in the first pass, but using more strin- 
gent criteria of 5 = 0.7 and SNR = 3.0. This ensures that 
fainter pixels around cosmic rays are also fiagged. 

In addition, the ACS /WFC exposures are significantly 
impacted by CTE degradation, which introduces sub- 
stantial deferred-charge trails extending along columns 
away from cosmic ray hits. Although this has already 
been accounted for to some e xtent b y means of the CTE 
algorithm described in Section 5.2.2[ some residual effects 
can remain in the images near the locations of bright 
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Table 4 

Astrometric Reference Catalogs Used for Each CANDELS field 



Field 


Area 


Telescope 


Filter 


Depth (5o- AB) 


Resolution 


Reference 


GOODS-N 


34' X 27' 


Subaru/Suprime-C am 


R 


26.6 


I'.'l 


Capak et al. (2004) 






HST/ACS 


F850LP 


27.4 


O'.'OB 


Giavalisco et al.(2004) 


GOODS-S 


34' X 33' 


ESO/2.2m WFI 


R 


25.5 


O'.'S 


Arnouts et al.(2001) 






HST/ACS 


F850LP 


27.4 


O'.'OB 


Giavalisco et al.(2004) 


COSMOS 


2° X 2° 


Subaru /Suprime-Cam 


i+ 


26.2 


0"9 


Capak et al. (2007) 






HST/ACS 


F814W 


27.2 


0^08 


Koekemoer et al. (2007) 


ECS 


1° X 1° 


CFHT/12k 


R 


24.7 


o'.'g 


Coil et al. (2004) 






HST/ACS 


F814W 


27.2 


o'.'os 


Davis et al. (2007) 


UDS 


1° X 1° 


UKIRT/WFCAM 


K 


24 


0'.'8 


Lawrence et al. (2007) 



cosmic rays. These are mitigated in MultiDrizzle by 
a subsequent rejection iteration whereby the cosmic ray 
masks are convolved with a linear kernel in the direction 
opposite to the read-out direction, and with a length of 
30 pixels, which is about the worst-case length of the 
CTE trails. The rejection is then carried out in this ad- 
ditional region, using the same parameters as were used 
in the second-pass iteration for pixels around the origi- 
nally flagged ones. 

For the WFC3/IR images, most of the cosmic rays 
are already rejected during the up-the-ramp sampling. 
However, there are occasional cosmic rays that are not 
fully removed, or warm pixels that are not accounted 
for in the dark file correction, so the WFC3 exposures 
for each filter are also passed through this step. For 
the ACS/WFC data, there are typically between two 
and four exposures per filter in a given orbit. Given 
the typical cosmic ray rate of ^ 1 — 2% during our ex- 
posure times, this means that for a 4-exposure depth, 
< 1 — 2 pixels can be expected to be hit by cosmic rays 
during all four exposure, while for two exposures this 
number increases to ^ 2000 — 6000 pixels would be af- 
fected by cosmic rays during both exposures. However, 
this is still only ~ 0.01 — 0.04%, meaning that ~ 1 out of 
every 100 small galaxies (~ 100 pixels in area) would be 
affected, losing ~ 1 — 4% of its pixels on average. 

5.5. Absolute Astrometric Calibration 

The absolute astrometric accuracy of any HST obser- 
vation is limited by the astrometric uncertainty of the 
primary guidestar, which for the GSC-2 is ^ 0'.'3 — O'.'S. 
In addition to a shift, these uncertainties also introduce 
errors in the knowledge of the orientation of the tele- 
scope, which is derived from the guidestar positions as 
well. All these uncertainties can be problematic for a pro- 
gram like CANDELS, where observations from many dif- 
ferent visits and epochs, taken with potentially different 
guidestars, need to be combined into a single set of im- 
ages, and where accuracy to levels better than ~ 0.1 pixel 
are demanded. We address this using a 2-stage tech- 
nique: firstly, obtaining the best possible astrometric so- 
lution for all the visits relative to one another by using 
a deep external ground-based catalog of the field which 
provides sufficient source density to match enough ob- 
jects on the HST images, and secondly by ensuring that 
the entire ground-based catalog / mosaic reference frame 
is registered to an absolute external astrometric reference 
coordinate system. 



For each of the five different CANDELS fields, appro- 
priate external astrometric reference catalogs are avail- 
able, which we list in Table |4] and describe here briefly. 
For GOODS-North, the astrometry is based on the 



Suba ru / Suprime-Cam R-band imaging (Capak et al. 
2004[ ), whose absolute astrometric reference frame has 
been registered to SDSS, 2MA SS, and the new deep 
VLA 20 cm survey of the field (Morriso n et al.||2010[ i. 
For GOODS-South, we use an R-band mosaic from the 
ESO-MPI 2.2m / Wide Field Imager (WFI) obtained 
as part of the ESQ Imaging Survey (EIS) of t he field 



(Arnouts et al. 



et al. 



2001[ ), registered to the GSC-2 ( [Lasker 

ToFtlie COSMOS field, the catalog used 

is based on CFHT/Megacam i* imaging, supplemented 
by deeper Subaru/ Suprime-Ca m /+ imaging (both de 



scribed in Capak et al. 2007), with absolute astrome 



try determined by regi stering to the VLA 20 c m survey 



of the COSMOS field ( Schinnerer et al. 
EGS field, we use the Die^2 CFHT/12' 



2004). For the 



ing of the field ( Coil et al. 



2004 



t mosaic i mag- 



Davis et al. 2007), ini 



tially registered to the USiNO-A2.0 catalog and subse- 
qucntly imp roved by matching to the VLA 20 cm survey 
of this field (Ivison et al. 2007). The absolute astromet- 



ric accuracy is generally ^ U".' 1 for each of these reference 
catalogs, limited primarily by residual uncertainties in 
the underlying reference frame. In addition, since exten- 
sive HST imaging is available for all these fields, which 
has already been registered onto the astrometric systems 
provided by these catalogs and are much deeper with 
better resolution, we use these HST data in determin- 
ing the astrometry for the new CANDELS data, and 
only use the ground-based catalogs in regions that are 
not covered by the existing HST data. Only the fifth 
field, UKIDSS/UDS, has insufficient existing HST cov- 
erage. Therefore our fundamental astrometric frame for 
this is t he UKIDSS K-band c atalog (Almaini et al., priv. 
comm.; Lawrence et al. 2007), which has been obtained 



by imaging the field with UKIRT/WFCAM a nd is reg- 



istered to t he VLA 20 cm survey of the field ( Simpson 



et al. 2006 ) to a similar level of astrometric accuracy as 
the other surveys . 

For all the exposures of a given instrument in each or- 
bit, an "astrometric detection image" is then produced 
by doing another run of MultiDrizzle to produce a sin- 
gle combined image (separately for each different instru- 
ment), for the tile covered during that orbit, containing 
all the exposures of all the filters for that instrument, and 
applying the cosmic ray masks that have been produced 
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Figure 5. Results showing the residual uncertainties in the positions of all the sources in the first epoch of GOODS-South as an example, 
where the plotted values indicate the difference in position for each source as measured from the new CANDELS data, compared with the 
position of the same source in the reference GOODS astrometric catalog, for the ACS/WFC images (left) and WFC3/IR images (right). For 
WFC3/IR the larger uncertainties are predominantly related to the broader PSF and larger pixel size, relative to ACS. The new GANDELS 
ACS data that overlap the previous GOODS- ACS data generally have smaller residuals than the WFC3/IR data. However, a subset of 
the new CANDELS AGS data lie outside the existing GOODS-AGS region, where we had to register the images to the ground-based WFI 
R-band data, and these sources have larger residuals (up to '^O'.'l per source). Given ~ 300 — 400 sources per tile, with the measurement 
errors on each source being drawn from this distribution, the overall accuracy on the absolute astrometric alignment of a given tile is 
therefore <0'.'005, once the uncertainties on the positions of all the individual objects have been combined in quadrature. Therefore, for 
both instruments, the resulting absolute astrometry for each orbit, after combining the astrometric information from all the sources, is 
accurate to better than < 0.1 pixel, sufficient to enable robust combination of data in overlapping regions as well as permitting cosmic ray 
rejection across multiple epochs once the images have been placed onto this common astrometric grid. 
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Figure 6. Bimicd vector residual plots showing the distribution of astrometric uncertainties across the GOODS-S field, once the images 
have all been registered onto the GOODS astrometric grid, again for the first epoch of GOODS-South, showing the ACS/WFC data (left) 
and WFG3/IR (right). The bimicd cells arc 40" in extent, where the size is chosen for display purposes to ensure a sufficient number of 
objects per cell while also providing a sufficient number of cells across the mosaic to show the general structure of the residuals. The two 
cameras arc offset from one another since the ACS/WFC data were obtained in parallel to the WFC3/IR data. Note also the sharp change 
in accuracy for the ACS data toward the south-west. This area lies outside the GOODS ACS v2.0 catalog, where we had to register the 
images to the ground-based WFI R.-band data. However, even in this region the accuracy remains reliable to <0'.'l for each individual 
source, which means that when an overall shift is computed for these exposures using the full sample of ^ 300 — 400 sources per exposure, 
the final shift is a<;curate to < O'.'OOS — O'.'Ol once the measurements from all the different sources are combined, thereby enabling robust 
alignment of ea/ch exposure relative to those from different orbits. 
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when combining the filters separately. This is motivated 
by the fact that the relative shifts for all the exposures, 
for all filters in a given orbit, have now already been 
determined to an accuracy of ^0.5— Imas. Therefore 
nothing further would be gained by now attempting to 
solve for shifts separately for different filters in a given 
orbit. Therefore, the exposures for all the filters are com- 
bined into a single image (separately for each different 
instrument), yielding one image for WFC3/IR, another 
for WFC3/UVIS, and another for ACS/WFC, for each 
orbit. A catalog is then produced from this single multi- 
filter image, which also has the advantage of providing 
increased depth and reducing the impact from cosmic 
rays compared with the individual filter images. 

All the sources in the multi-filter catalog for each orbit 
are then matched to the sources in the relevant portion of 
the external catalog, using a number of iterative steps. 
The first iteration uses a relatively large tolerance (up 
to a few arcseconds) and only the brightest ^20 — 30 
sources in each image, in order to determine the domi- 
nant terms in the shifts for right ascension and declina- 
tion. Once these have been accounted for, several addi- 
tional iterations are carried out using the full catalog of 
sources in each image, using progressively tighter match- 
ing tolerances down to C/.'l and solving for the residual 
remaining shifts as well as the rotation errors due to 
the uncertainties in guidestar position. For all visits, 
^ 300 — 400 sources are typically then matched at the 
faintest levels and tightest tolerances between the HST 
MultiDrizzle-combined image and the reference cata- 
log. 

The results of this procedure are shown in Figures [5] 
andjel for the ACS/WFC and WFC3/IR images of the 
GOODS-S field. Figure [5] shows the distribut ion in right 
ascension and declination offsets between the reference 
catalog positions of the sources and those measured on 
the new CANDELS data, after having solved for the as- 
trometry as described above, for all sources in the first 
epoch of the GOODS-South field by way of example. In 
general, the offset between the new measured position 
and the catalog position is < O'.'l per object, where these 
residual differences are due to differences in sensitivity 
between the new and the old data, as well as morphologi- 
cal differences due to different filters, and astrometric un- 
certainties in the reference catalog. Since the ACS /WFC 
and WFC3 /IR images for each orbit typically contain at 
least ^ 300 — 400 sources, each with a residual < O'.'l, this 
means that the overall absolute astrometric solutions for 
the exposures in each orbit are known to better than 
<0'.'005 (i.e., ~ 0.1 pixel for ACS), once the uncertain- 
ties on the positions of all the individual objects have 
been combined in quadrature. 

Figure [6] is based on the same data as in Figure [5] but 
this time showing the residuals as a function of position 
across the field, where each vector indicates the mean 
residual in a grid of cells, each 40" on a side, where the 
size is chosen for display purposes to ensure a sufhcicnt 
number of objects per cell while also providing a sufhcient 
number of cells across the mosaic to show the general 
structure of the residuals. The residuals are generally be- 
low ~],0'.'05 for most sources across the GOODS-South 
region, except for the CANDELS fields in the south- 
western corner which lie outside the previous GOODS- 
ACS coverage. For those fields, the only available refer- 



ence catalog positions are those from the ground-based 
2.2m ESO WFI R-band catalog previously discussed, and 
therefore the residuals per object are somewhat larger in 
that region. However, even for those tiles, the overall as- 
trometric solutions for each orbit better than ^0.1 pixel 
for ACS once all the objects in all the exposures are taken 
into account, thereby enabling these exposures to be ac- 
curately aligned with the rest of the data. 

We are continuing to investigate further improvements 
to these technique as we accumulate more data on the 
CANDELS fields, including the possibility of solving si- 
multaneously for the positions of all sources, on all expo- 
sures, for all filters and all pointings, using the ground- 
based catalog as an external constraint. This may have 
the additional advantage of solving for the small resid- 
ual scale changes from one exposure to the next, as well 
as any possible changes in rotation during an orbit, and 
perhaps even long-term changes in the skew terms of the 
detector which might not yet be captured in the distor- 
tion models. However, the current approach yields an 
absolute astrometric accuracy on the order of ~ 5 mas, 
or ^0.1 pixel for ACS/WFC, between images obtained 
with different guidestars. This is sufficient to prevent reg- 
istration problems between different HST observations 
of the same object in overlapping images from different 
visits or epochs and thereby preserve the morphological 
properties of galaxies in the survey. 

5.6. Satellite Trails, Optical Ghosts, and Crosstalk 

In addition to cosmic rays, several other artifacts are 
also present in the images and are best identified vi- 
sually and masked by hand. These include trails from 
bright satellites, optical filter ghosts from bright stars, 
and anomalous persistence signals that might not have 
been identified in other ways. To this end, a data qual- 
ity validation team has been assembled from across the 
CANDELS collaboration, who access all the individual 
exposures and MultiDrizzle-combined images by means 
of a set of webpages that enable them to visually inspect 
graphical displays of the images, as well as to retrieve 
the images and submit data quality input results that 
are subsequently collected and used to identify and cor- 
rect anomalous images. 

Due to the relatively low altitude of the orbit of HST, 
several percent of exposures are affected by the passage of 
satellites across the field of view during the exposure. In 
addition, the WFC3/IR detector can exhibit persistence 
at the location of the satellite trail for up to several ad- 
ditional orbits if the satellite is bright enough. All expo- 
sures that are affected by this are identified visually, and 
a software script is then run to mask the satellite trail 
in the undistorted image and subsequently transform the 
mask back to the distorted frame of the input exposure. 
These masks are then included along with the cosmic 
ray masks in identifying pixels to be excluded from the 
combined mosaics for each individual epoch. 

Optical ghosts from stars are also visually identified 
and masked on the affected exposures. In this case the 
masks are not used in creating the combined image for 
each individual epoch, since the dither pattern used is 
sufficiently small that the stellar ghosts do not move by 
much, which would then result in large holes in the result- 
ing image. Instead, these masks are used when combin- 
ing multiple epochs, since the locations of optical ghosts 
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in one epoch are generally covered by pixels from another 
epoch that are not affected by ghosts. 

Finally, an electronic effect that is present particularly 
in CCDs, and to a lesser extent in the WFC3/IR detec- 
tor, is electronic crosstalk from bright sources, producing 
a region of somewhat lower flux that is located symmet- 
rically in another quadrant. For ACS/WFC this effect is 
substantially mitigated by our use of GAIN=2 for all ex- 
posures and is not present at any significant level in the 
CANDELS data. Specifically, the inter-CCD cross-talk 
is negligible at GAIN=2, while the intra-CCD cross-talk 
of ^ 0.007% is corrected by the current version of the bias 
de-striping code. For WFC3/UV1S, however, it is quite 
noticeable particularly for stars that are bright enough 
to have a charge-bleeding column of bright pixels, lead- 
ing to a corresponding low-valued section of columns on 
the opposite quadrant. The brightest galaxies in these 
fields also introduce this effect, with the affected region 
being more diffuse in that case. For WFC3/IR the eifect 
is observed for only the very brightest stars and consists 
simply of a region with a relatively minor decrease in 
fiux corresponding approximately to the area of satura- 
tion of the star. In our CANDELS imaging pipelines, we 
are able to flag the worst of the WFC3/UVIS crosstalk 
by identifying bleeding charge-trail columns in bright 
sources and then replicating a pixel mask for the corre- 
sponding area on the opposite side of the detector. These 
pixels are then excluded when creating the combined im- 
ages for each individual epoch, since the dither offsets are 
large enough to ensure that such regions are generally 
sufficiently covered by other pixels that are not affected 
by this issue. 

5.7. Residual Background Subtraction 

Prior to combining all the exposures into a flnal mosaic, 
it is necessary to remove all background emission due 
to non-astronomical sources, for example the sky back- 
ground from earthshine, which is relevant to all three 
cameras, as well as zodiacal light and low-level thermal 
background emission, which are more applicable to the 
WFC3/IR camera. Because these vary with time, their 
contribution to the pixel counts would lead to photo- 
metric errors in the final count-rates of sources if their 
relative differences between one exposure and the next 
are not removed. Moreover, the early versions of the 
WFC3/IR flatfields contained significant low-level large- 
scale residuals, ~3 — 4%, and we first embarked on an 
effort to improve these by constructing skyflats based on 
empty regions in all the CANDELS data that we had to 
date. However, new skyflats have now been released by 
the WFC3 team that provide improved flatfield correc- 
tions to ~1% for WFC3/IR, and after we incorporated 
these into the CANDELS products we found that they 
significantly reduced the large-scale structure across the 
images, as well as improving the consistency of photo- 
metric calibration across the detector. 

In addition, some of the most severe instances of sky 
background emission can come from cases where the tele- 
scope is observing close to the limb of the bright earth, 
with backgroimds increasing significantly for bright earth 
limb angles below ~ 35 — 40° , and in some cases being 
non-uniform across the detector. For the WFC3/UVIS 
and ACS/WFC exposures not much can be done to mit- 
igate this in any given exposure, since each CCD expo- 



sure is just a single integration. For WFC3/IR, however, 
we have been investigating the possibility of excluding 
specific reads near the beginning or end of the MULTI- 
ACCUM sequence, since the reads are sampled at either 
50s or IGOs and often only a few of them arc affected 
by bright earth limb emission near the beginning or end 
of the orbit, while the remainder of the reads are ob- 
tained at higher limb angles with a more nominal back- 
ground. Thus excluding the high-background ones can 
potentially improve the S/N (since the loss in S/N from 
excluding some reads can be compensated for by the in- 
creased S /N due to a lower background in the remaining 
ones). The exposures obtained so far generally do not 
show a sufficient degree of improvement to warrant this, 
and therefore we are currently retaining all the reads, but 
this technique remains available as an option for any fu- 
ture observations that may potentially be significantly af- 
fected by unusually bright sky backgrounds at low earth 
limb angles. 

For all exposures, once all the low-level residual struc- 
ture has been removed, the mean background level 
is then determined from a masked version of the ex- 
posure, which is constructed as follows. Firstly, the 
clean MultiDrizzle-combined image for each visit is run 
through a source-detection step to create a mask of all 
sources in the images, which is done by smoothing the 
image and applying a sigma-clipping threshold for pixels 
that exceed the mean background level. This is neces- 
sary because many sources display faint isophotes that 
reach into the noise of the image and are not excluded 
when doing a simple sigma-clip on the image, thereby 
contributing to the pixel statistics at faint levels and bi- 
asing the resulting sky estimate. Smoothing the image 
effectively increases the significance of these faint outer 
isophotes and thereby enables these pixels to be flagged 
to much fainter levels than in the unsmoothed image, 
thereby reducing their impact on the final sky estimates. 
Typically a total of ~ 20 — 30% of all the pixels in the 
images are excluded in this way, leading to a significant 
reduction in the impact of these pixels on the final pixel 
statistics and a much more stable background sky esti- 
mate. 

The statistics of the remaining pixels are still slightly 
biased toward positive values, which is simply a conse- 
quence of the much larger number of faint sources be- 
low the detection thresholds of the images, whose com- 
bined signal adds a slight net positive bias to the pixels. 
To these are added other faint positive signals such as 
faint tails from cosmic rays and deferred-charge trails 
from CTE which in both cases may be below the detec- 
tion thresholds for the algorithms that are responsible 
for identifying and masking these. Ultimately, however, 
these residual effects are approximately uniform from one 
exposure to the next and represent a net DC component, 
and all that changes between exposures is the varying 
sky, zodiacal and thermal background. Therefore, the 
background level on all the remaining unmasked pixels 
on each image (~ 70 — 80% of the total number of pixels) 
is determined by doing an iterative sigma-clipped fit to 
these pixels and determining its mean, which is then used 
as the final background offset value and subtracted from 
each individual exposure to yield an image containing 
only the counts from the sources. 
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Table 5 

World Coordinate System Information for each CANDELS HST Mosaic 



Field 


Mosaic Tangent Point 


Instrument / Camera 
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5.8. Final MultiDrizzle Mosaic Combination 
5.8.1. Inverse Variance Images 

In preparation for the final step of combining all 
the exposures for each filter into a single mosaic using 
MultiDrizzle, our pipelines first create for each expo- 
sure a corresponding inverse variance image, which con- 
tains all the "intrinsic" error terms associated with each 
pixel (including noise from accumulated dark current, de- 
tector read-out, and photon noise from the background 
as modulated multiplicatively by the flatfield and the de- 
tector gain), but not the additional "extrinsic" Poisson 
noise from the astronomical sources in the image. T his 
approach was described by Casertano et al. (2000) for 
the Hubble Deep Field South and has also be en used 
in GOODS (Giava lisco eral|[2004| , COSMOS ( |Scoville 
et al.|2007| ,Koekemoer et al.|2007[ f, AEGIS (D avis et al. 
2007; Mewman et al. 2011, in prep.), UDF (IBeckwith 



et al. 20061, and other projects, and i s implemented as 
a rou tine option in MultiDrizzle (see Koekemoer et al. 
2002 1. We provide here a brief description ot it as ap- 
plied to the new instruments that are used to observe 
the CANDELS fields, but refer the reader to the above 
references for further details. 

When combining images where each pixel value has 
noise associated with it, the optimal approach involves 
weighting each pixel according to the inverse square of 
its noise. In particular, it can be demonstrated that the 
appropriate quantity to use consists of the background 
noise only; it should not include the noise associated with 
the emission from the object since that would lead to a bi- 
ased estimation of the true flux in the pixel and result in 
photometric errors. Our CANDELS pipeline constructs 
inverse variance images for all the exposures obtained 
with WFC3/IR, WFC3/UVIS and ACS/WFC, using es- 
sentially the same formalism as described in Casertano 
et al. (2000) for the original WFPC2 Hubble Deep Field 
South observations, but modified to take into account the 
differences in the data formats for the ACS and WFC3 
detectors. In particular, these detectors are calibrated 
to electrons per second and have been corrected for the 
detector gain, while the WFPC2 data are in DN and 
therefore still need to have the gain included in the cal- 
culation. Specifically, the WFC3/1R images are in elec- 
trons per second while WFC3/UVIS and ACS/WFC are 
in units of total electrons; however, the input inverse vari- 



ance images for all of them are calculated to be in units 
of counts per second, since the ACS images are scaled 
to be in electrons /second while drizzling and since that 
is also what the output units of the mosaics are in all 
cases. The formula that is used to calculate the inverse 
variance images, in units of (e/s)~2, is given by: 



Inverse Var. = 



iftf 



(D + fB) + a 



read 



where / is the inverse flatfield (as defined in the conven- 
tions for the flatfield reference files used in calibration), t 
is the exposure time (in seconds), D is the total accumu- 
lated dark current signal during the exposure, B is the 
total accumulated background level during the exposure, 
and (Tread IS the read-out noise, with all three of the latter 
quantities being in units of electrons. 

5.8.2. Pixel Scale and Final Drizzle Parameters 

The pixel scale for the output mosaics is driven by 
the detector plate scale and pixel size, together with the 
full width half max (FWHM) of the point spread func- 
tion (PSF) produced by the telescope optics. At the 
wavelengths of the WFC3/IR F105W to F160W observa- 
tions, the HST PSF has a FWHM - 0'.'12 - O'.'IS, which 
is subsequently convolved by the 0'.'128 WFC3/IR de- 
tector pixel scale. Hence, the best PSF that could be 
recovered (without deconvolution) , even in the ideal sce- 
nario of combining images using interlacing, which would 
minimize additional convolutions, still has a FWHM 
'-^ 0'.'17 — 0'.'19 in the final images. We choose an output 
pixel scale of O'.'06/pixel for the final WFC3/IR mosaics, 
providing adequate sampling of the PSF. 

Having determined an output pixel scale facilitates the 
choice of the other relevant Drizzle parameter, namely 
pixf rac, which defines how much the input pixels are re- 
duced in linear size before being mapped onto the output 
grid (see the basic descr iption of Drizzle parameters in 
Koekemoer et al. 2002 1. The initial PSF produced by 
the HST optics is first convolved by the detector pixel 
size when imaging the sky, then a second time by the 
rescaled detector pixel size when mapping onto the out- 
put grid (thus applying the pixf rac parameter), then 
is convolved a final time by the output pixel scale. The 
impact of the second convolution can thus be minimized 
by setting pixf rac to a sufficiently small value; setting 
it to would remove this convolution, corresponding to 
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Figure 7. Photometric comparison of magnitudes measured from the CANDELS WFC3/IR imaging of the UKIDSS/UDS field, relative 
to the photometry of the same objects from the publicly released UKIDSS J and H data. Since the filter bandpasscs are somewhat different , 
the UDS photometry has been transformed to the WFC3/IR system using stellar spectral libraries (Pickles,1998, ,Gunn fc Stryker|l"983l l, 
thereby enabling a direct comparison of any remaining possible offsets. The agreement is generally robust to better tnan a tew percent, in 
agreement with the current published photometric performance of WFC3. 



pure interlacing (as was done for the UDF as described 
in Beckwith et al. 2006, where up to 144 exposures were 
available) , but for a program like CANDELS with a lim- 
ited number of exposures, we set pixfrac = 0.8 in order 
to avoid introducing too much variation between pixels 
in the corresponding weight images. 

For the ACS/WFC and WFC3/UVIS exposures, the 
PSF is considerably sharper due to both the shorter 
wavelengths and the smaller pixel scale on the detector; 
thus we choose an output mosaic pixel scale of O'.'03/pixcl 
for these (which has become the standard for GOODS, 
UDF, and most other large HST ACS surveys), and also 
set pixf rac — 0.8 based on similar considerations as with 
the WFC3/IR data. The resulting PSF, after taking into 
account the relevant convolutions, is ~ O'.'O? — O'.'ll in the 
final images, across the range of UV/optical filters that 
we are using and is very well sampled by our chosen 
O'.'03/pixel scale. 

The final pass of MultiDrizzle is then run using the 
above output pixel scale and pixfrac settings and ap- 
plying the inverse variance weight image associated with 
each exposure, which contains the full set of masks from 
cosmic rays, bad pixels, satellite trails, and other blem- 



ishes in the detector. In each case the images for all filters 
are drizzled onto a common tangent plane projection on 
the sky, which is defined to match the existing ones where 
known, to facilitate a direct comparison with pre-existing 
data on these fields. Four of the fields (GOODS-North, 
GOODS-South, COSMOS, and ECS) all have an exist- 
ing tangent plane point already defined, which we adopt 
for the CANDELS mosaics as well. For the UDS field, we 
adopt a common tangent plane point designed to satisfy 
the current surveys on that field. The World Coordinate 
System (WCS) properties of all the CANDELS mosaics 
are shown in Table [S) 

Figures [T0| to [25| show the final combined mosaics ob- 
tained to date for the GOODS-S and UDS fields, in- 
cluding the WFC3/IR and UVIS observations, which are 
all obtained as primary exposures. Also shown are the 
ACS/WFC exposures, which are obtained in parallel and 
are therefore offset from the WFC3 pointings. In each 
case we show a color mosaic representing the various fil- 
ters obtained in each given instrument (either WFC3/IR 
or WFC3/UVIS-HACS), together with the inverse vari- 
ance weight images which show the total extent of the 
coverage and degree of overlap between the WFC3/IR 
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Figure 8. Plot of measured FWHM values for all the sources in the CANDELS UDS field, for all four filters (ACS/WFC F606W and 
F814W, and WFC3/IR F125W and F160W). Stars that have been manually inspected and verified to be valid, isolated point sources with 
no significant defects are identified in red. The mean FWHM for these stars is shown as a solid line, with dashed lines indicatring the 
1-sigma standard deviation boundaries in each case. 
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Figure 9. Plot of measured FWHM values for the subset of manually verified unresolved stars in CANDELS UDS field, for all four filters 
(ACS/WFC F606W and F814W, and WFC3/IR F125W and F160W), as a function of position across the field. The plots indicate the 
mean FWHM with a solid line and the standard deviation boundaries using dashed lines. Generally the stars are well behaved across the 
field, with no significant trends evident as a function of field position. We note that the somewhat higher scatter of the WFC3/IR points is 
due primarily to the increased importance of pixellation due to its larger detector pixel size, thus the ACS measurements are our primary 
diagnostic and these do not reveal any significant trends as a function of position across the field. 
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pointings as well as the WFC3/UVIS and ACS/WFC im- 
ages. The WFC3/IR exposures successfully adjoin one 
another with minimal overlap while maximizing the area, 
as designed. However, the larger ACS/WFC exposures 
overlap each other significantly, given that the offsets 
from one pointing to another are governed by the size 
of the WFC3/IR detector. The ACS/WFC exposures 
are therefore effectively doubled at most locations across 
the field. 

Three months after the observations of a given CAN- 
DELS epoch are completed, we release the calibrated mo- 
saics to the public via the STScI archivep^ including the 
drizzled science mosaics as well as the inverse variance 
weight files that describe the noise associated with each 
pixel. Updates on the data obtained will be provided at 
the primary CANDELS project website]^ as the survey 
progresses. Subsequent papers will present details on the 
catalogs and other measurements derived from the CAN- 
DELS data. 

5.9. Photometric Validation 

We have carried out a series of photometric tests be- 
tween the HST and ground-based imaging in fields for 
which both types of datasets exist, in order to quan- 
tify the level of agreement between photometry derived 
from the two types of data, with emphasis on the new 
WFC3/IR imaging. In particular we examined the 
UKIDSS/UDS field since this is the first CANDELS 
field to have been completed, and extensive ground-based 
imaging exists in the near-IR which we can compare with 
our new data. The total WFC3/IR system throughput in 
the F125W and F160W filters, including the filter trans- 
mission curves, the detector sensitivity, and the optical 
telescope assembly response, is somewhat different from 
the ground-based system, therefore we carried out a rig- 
orous conversion between the two bandpasses using the 
following procedure. 

The first step involved identifying all the unresolved 
objects in the UDS field that were covered by CANDELS 
F125W and F160W as well as UKIDSS J and H, also en- 
suring that the objects were unsaturated in all cases and 
in a magnitude range that is well covered by the UDS 
data. The UKIDSS magnitudes were then converted to 
the WF C3/IR system using a set of stellar spectral li- 
braries ( |Pickles|[l998{ |Gunn fc Stryker||1983"| to obtain 
the required filter conversions, using detailed existing in- 
formation about the filters and total system throughput 
in each case, following a methodology similar to previ- 
ous multi-band catalog studies that have involv ed obser 
vations from several different instr uments (e.g., Grazian 
et al.||2006[ [Castellano et al.| |2010[ ). We plot the result- 
mg comparison in I^'igure [7j I'he agreement is generally 
tight and shows no significant systematics, with the ze- 
ropoints agreeing to better than a few percent, which is 
also in good agreement with what has been reported for 
the performance of WFC3 ( [MacKenty eTa l. 2010 Q We 
have also verified that initial checks on the data obtained 
so far for the other partially completed CANDELS fields 
yield similar results. 



^® http: / / archive.stsci.edu / prepds / candels / 
http://candels.ucolick.org/ 



See also http:/ /www. stsci.edu/hst/wfc3 



5.10. Point-Spread Function Validation 

The quality of the PSF is important to verify, both 
in terms of the scientific utility of the data (including 
morphological measurements and signal-to-noise limits 
for compact or unresolved sources) as well as the quality 
of the data (providing a good diagnostic of the alignment 
accuracy between exposures and the cosmic ray rejection 
quality for compact and unresolved sources). 

Within the CANDELS team, we have compiled an 
extensive list of unresolved sources in each pointing 
obtained to date (verified by manual inspection), and 
we use these here to diagnose the PSF quality of the 
UKIDSS/UDS dataset, since this is the first field to 
be completed and has the most extensive CANDELS 
WFC3/IR and ACS/WFC imaging. In Figure M we 
show the measured PSF of all the sources in this neld, 
for all four filters (ACS/WFC F606W and F814W, and 
WFC3/IR F125W and F160W), and identifying in red 
the stars that have been manually verified as being valid 
point sources that are not saturated and are clearly iso- 
lated. 

Generally the unresolved sources are very well behaved, 
until reaching ~ 18th magnitude where saturation be- 
comes an issue. We note that the ACS FWHM are 
slightly broader than might be expected, because all 
data in this particular dataset have been drizzled to a 
O'.'06/pixel scale, primarily to facilitate cross-image cat- 
aloging and related work. However we have verified that 
images drizzled to a finer scale (O'.'03/pixel) recover the 
- 0'.'08 - 0'.'09 FWHM generally achieved for ACS imag- 
ing. The WFC3/IR data display a broader FWHM pri- 
marily as a result of the broader HST PSF in the IR 
as well as the larger ^0'.'13 detector pixel scale, which 
results in more dither-related pixellation due to the rel- 
atively small number of dither offset positions used. 

We also examine the measured PSF for this subset of 
unresolved sources as a function of position across the 
field, and show the result in Figure [9) This test is useful 
for identifying regions of the image where misalignment 
or cosmic ray rejection problems might cause a net up- 
ward shift of the FWHM by a significant amount. How- 
ever, it is evident that across the field the PSF values re- 
main well behaved, especially for ACS which is more sen- 
sitive to CR rejection issues resulting from misalignment, 
and that no significant systematics are present across the 
field. 

6. SUMMARY 

We have described the i/ST observational imaging data 
products and processing pipelines for the Cosmic As- 
sembly Near-infrared Deep Extragalactic Legacy Survey 
(CANDELS), a 902-orbit F^TMuffi-Cycle Treasury pro- 
gram aimed at documenting the first third of galactic 
evolution from redshift z ^ 8 to 1.5 via deep imaging of 
more than 250,000 galaxies with the Hubble Space Tele- 
scope, together with new constraints on the use of SNe la 
as tracers of dark energy to higher redshifts than previ- 
ously studied. The survey covers five well-studied extra- 
galactic fields, namely GOODS-North, GOODS-South, 
COSMOS, EGS and UKIDSS/UDS, targeting each one 
predominantly with WFC3/IR as prime and ACS/WFC 
in parallel, together with some WFC3/UVIS observa- 
tions where necessary for the science goals. The data are 
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all recalibrated and combined using our own pipelines as 
described here. Three months after the observations of 
a given CANDELS epoch are completed, wo release the 
calibrated mosaics to the public via the STScI archive. 
At the time of writing, we have already released the first 
few GOODS-South epochs, and the full UDS campaign. 
We strongly encourage the astronomical community to 
make use of the CANDELS data to advance their own 
research. 
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APPENDIX 

ADDITIONAL HST DATA IN THE CANDELS FIELDS 

In addition to the new CANDELS data described in this paper, several of these survey fields also contain significant 
previous _H5T investment, in many cases with instrument/filter combinations that provide further scientific leverage 
when combined with the CANDELS data. These surveys include GOODS (Giavahsco et al. 2004), GEMS (Rix et al. 
2004), UDF (Beckwith et al. 2006; Oesch et al, 2007; Bouwens et al. 2010), WFC3-ERS2 (Windhorst et al. 2011), 
COSMOS (Scoville et al. 2007), and ECS/AEGIS (Davis et al. 2007; Newman et al. 201f), as weU as a variety of 
other programs that have obtained data in various locations on these fields. We note that the HST datasets for these 
programs are all already public, in many cases with high-level science products already released, although in most cases 
the calibrations have improved since the data products were originally released. Thus we are reprocessing these data, 
combining them with newer data using the same CANDELS filters where appropriate and creating summed mosaic 
images that have been uniformly processed using updated calibration files and CANDELS methods as described in 
this paper. We note also that the exposure maps and sensitivity estimates of Grogin et al. (2011) take these prior HST 
programs into account. These datasets are summarized in Table [6j showing the existing (non-CANDELS) broadband 
filter coverage in ACS/WFC and WFC3/IR and UVIS, as well as the relevant references where available. 

The datasets are all retrieved from the archive and are recalibrated using the current generation of reference files, 
including all the new fiatfield, dark current, bias and geometric distortion files, as well as improved photometric 
calibrations. We note that some of the earlier ACS data do not suffer from the same electronic issues as the post-SM4 
ACS data, most notably the bias striping effect, and CTE degradation is also much less for the first year or two of ACS 
operations (since CTE degradation gradually worsens over time). However, the flatfields, darks, bias and distortion 
reference files are generally much improved, both for the older ACS data and for the newer WFC3 data, and are 
therefore applied to these reprocessed datasets. These reprocessed datasets will be delivered and served out along with 
the final combined CANDELS datasets once all the data have been obtained for the survey. 



Table 6 

Additional ff5T ACS/WFC and WFC3/IR,UVIS Broad-band Filter Coverage in the CANDELS Fields 



Field 


ACS/WFC Filters 


WFC3/UVIS Filters 


WFC3/IR Filters 


GOODS-Ni 


F435W F475W F606W F775W F814W F850LP 




F140W 


GOODS-S^ 


F435W F475W F606W F775W F814W F850LP 


F225W F275W F336W 


F098M F105W F125W F140W F160W 


COSMOS^* 


F814W 


F300W 


F140W F160W 


EGS* 


F606W F814W 




F140W 


UDS 









1 Giavalisco et al. (2004); Malhotra et al. (2005); Rhodes et al. (2005); Riess et al. (2004) 

2 Giavalisco et al. (2004); Beckwith et al. (2006); Bouwens et al (2010); Malhotra et al. (2005); Rhoads et al. (2005); Riess et al. 
(2004); Rix et al. (2004); Oesch et al. (2007); Thompson et al. (2005); Windhorst et al. (2011) 

3 Scoville et al. (2007); Franx et al. (2008); Scarlata et al. (2007); van Dokkum et al. (2006) 

" Davis et al. (2007); Newman et al. (2011, in prep.); Franx et al. (2008); van Dokkum et al. (2006) 
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03^ 3,3'" CC^ J'::' JO" Ib^ Cu' y'"' 4':: 

Right Ascension (J2000) 



GOODS-S-Deep Epoch 1 




03^ 33"^ 00^ 32"^ 45^ 30^ 15^ 00^ Si"" 45^ 

Right Ascension (J2000) 

Figure 10. Images showing the prime CANDELS WFC3/IR dataset for the first epoch obtained on the GOODS-S field (GOODS-S-Deep 
Epoeh 1). This dataset also includes the first CANDELS test orbit (obtained in April 2010), to the north-west. The top panel shows a color 
composite of the WFC3/IR F125W and F160W images after combined mosaics were created for each filter separately using MultiDrizzle, 
while the bottom panel shows the corresponding weight images, which are in units of inverse variance. The F125W data are shown in blue 
and the F160W data are shown in red. Regions containing bad pixels (such as the circular "death star" region) are set to and thus have no 
weight in this single epoch dataset, in which the dither offsets were not yet large enough to move over such features. The overlap between 
pointings was chosen to be just large enough to provide contiguous coverage while also maximizing total area covered. Occasional tiles are 
intentionally tilted or offset to enable appropriate guide stars to be selected. The rectangular outlines indicate the nominal boundaries of 
the existing GOODS coverage. 
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GOODS-S-Deep Epoch 1 
WFC3/UVIS F350LP, ACS F814W 




03 3,3 CC 



Right Ascension (J2000) 



GOODS-S-Deep Epoch 1 
WFC3/UVIS F350LP, ACS-FSMW 




-28° 00 



03 33 00^ 32 45^ 30' 15' 00' 31 45' 

Right Ascension (J2000) 

Figure 11. Images showing the prime CANDELS WFC3/UVIS and parallel ACS/WFC datasets for the first epoch obtained on the 
GOODS-S field (GOODS-S-Decp Epoch 1). This dataset also includes the first CANDELS test orbit (obtained in April 2010), to the 
north-west. The top panel shows a color composite of the WFC3/UVIS F350LP and ACS/WFC F814W images after combined mosaics 
were created for each filter separately using MultiDrizzle, while the bottom panel shows the corresponding weight images, which are in 
units of inverse variance. The F35DLP data are shown in blue and the F814W data are shown in red. The WFC3/UVIS data consist of 
only one exposure per pointing and therefore still contain a large number of cosmic rays and other defects. Regions containing bad pixels 
are masked where necessary. Note also that, in general, the overlap between ACS pointings is sufficient to provide approximately twice the 
depth of a single pointing across much of the ACS area. The rectangular outlines indicate the nominal boundaries of the existing GOODS 
coverage. 
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GOODS-S-Deep Epoch 2 

WFC3/IR n25w, neow 




03^ 3,3'" CC^ J'::' JO" Ib^ Cu' y'"' 4':: 

Right Ascension (J2000) 



GOODS-S-Deep Epoch 2 




03^ 33"^ 00^ 32"^ 45^ 30^ 15^ 00^ 31"^ 45^ 

Right Ascension (J2000) 

Figure 12. As for Figure [To] but for GOODS-S-Deep Epoch 2. 
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GOODS-S-Deep Epoch 2 
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Figure 13. As for Figure [TT] but for GOODS-S-Deep Epoch 2. 
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GOODS-S-Deep Epoch 3 
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Figure 15. As for Figure [TT] but for GOODS-S-Deep Epoch 3. 
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GOODS-S-Wlde Epoch 1 
WFC3//UVIS F350LP, ACS F814W 
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Figure 17. As for Figure [TT] but for GOODS-S-Widc Epocii 1. 
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GOODS-S-Deep+Wlde Epochs 1-4 + ERS2 




03^ 3,3'" CG^ A'::' jif Ib^ Cu' 3''" 4':: 
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Figure 18. As for Figure [To] but tiiis time sliowing tlie full accumulated CANDELS dataset on the GOODS-S field so far (including 
GOODS-S-Deep Epochs 1, 27^ as well as GOODS-S-Wide Epoch 1), together also with the WFC3/IR ERS2 F125W and F160W dataset, 
since the CANDELS observations were designed to abut this dataset. 
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03^ 3,3'" CC^ A'::' jif' Ib^ Cu' 3''" 4':: 

night Ascension (J2000) 



GOODS-S-Deep+Wide Epochs 1-4 
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Right Ascension (J2000) 

Figure 19. As for Figure [IT] but thiis time sliowing tlie full accumulated CANDELS dataset on the GOODS-S field so far (including 
GOODS-S-Deep Epochs 1, 27^ as well as GOODS-S-Wide Epoch 1). 
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Figure 20. Images showing the prime CANDELS WFC3/IR dataset for the first epoch obtained on the UKIDSS/UDS field (UDS Epoch 
f). The top panel shows a color composite of the WFC3/IR F125W and F160W images after combined mosaics were created for each 
filter separately using MultiDrizzle, while the bottom panel shows the corresponding weight images, which are in units of inverse variance. 
The F125W data are shown in blue and the F160W data are shown in red. Regions containing bad pixels (such as the circular "death 
star" region) are set to and thus have no weight in this single epoch dataset, in which the dither offsets were not yet large enough to 
move over such features. Other masked features are largely satellite trails. The overlap between pointings was chosen to be just large 
enough to provide contiguous coverage while also maximizing total area covered. Occasional tiles are intentionally tilted or offset to enable 
appropriate guide stars to be selected. Exposures are shorter at the east end of the mosaic in order to accommodate F350LP exposures 
during the same orbit. 
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Figure 21. Images showing the parallel CANDELS ACS/WFC dataset for the first epoch obtained on the UKIDSS/UDS field (UDS 
Epoch 1). The top panel shows a color composite of the ACS/WFC F606W and F814W images after combined mosaics were created for 
each filter separately using MultiDrizzle, while the bottom panel shows the corresponding weight images, which are in units of inverse 
variance. The F606W data are shown in blue and the F814W data are shown in red. Regions containing bad pixels (such as satellite trails) 
are masked where necessary. Note also that, in general, the overlap between ACS pointings is sufficient to provide approximately twice the 
depth of a single pointing across much of the ACS area. 



36 



Koekemoer et al. 




Figure 22. As for Figure [20] but for UDS Epoch 2. 
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Figure 23. As for Figure [21] but for UDS Epoch 2. 
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Right Ascension (J2000) 



Figure 24. As for Figure [To| but sliowing thie full, final accumulated CANDELS dataset on tlie UDS field (including epochs 1 and 2). 
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Right Ascension (J2000) 



Figure 25. As for Figure [TTj but sliowing ttie full, final accumulated CANDELS dataset on tlie UDS field (including epochs 1 and 2). 



